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Der Rechenbereich einer numerischen Rechenmaschine kann dadurch auf Funktionen ausgedehnt werden, 

- daß die Maschine mit einem oder mehreren Speicherwerken versehen wird, in denen Funktionswerte für 

einzelne grob abgestufte Argumente und Interpolationsfaktoren gespeichert sind, aus welchen die Funktions- 

werte für beliebige Argumente durch ein- oder mehrstufige maschinelle Interpolation gebildet werden. Die 

technischen Lösungsmöglichkeiten werden angedeutet und die erreichbare Genauigkeit untersucht. Schließ- 
lich wird das Interpolationsverfahren mit anderen Verfahren verglichen. 

The range of action of a numerical compulalion machine can be ewpanded to functions by furnishing 
the machine with one or several storage mechanisms storing the values of the functions for large intervals 

. of the argument, together with interpolation factors. The values for any argument are obtained by mechanical 
interpolation by one or several sieps. T'he technical possibilities of construction are hinted, and the accuracy 
obtainable is ewamined. Finally, the method of interpolation is compared with other methods. 

Le domaine d’efficacit& d’une machine & calcul numerique peut Etre dtendu a fonctions, si on munit 
la machine d’un ou de plusieurs mecanismes capables d’enregistrer les valeurs de fonetions ü grands inter- 
valles de l’argument et des facteurs a’interpolation; on obtient les valeurs des fonctions pour un argument 
arbitraire par interpolation machinale a un seul ou ü plusieurs pas. Allusion est faite aux possibilites 
techniques de construction, et V’exactitude, qu’on peut atteindre, est ewaminee, Enfin, la methode d’inter- 
polation est comparee a d’autres methodes. 

O6dJacTb IIPHMeHeHHUA YUCHOBOH CYeTHOÜ MAIIMHBI MoHcer ÖBITB PacıpocrpaHena Ha 
BbIYMIICeHHe PYHEIMÜ TPM MOMOIM OMHOTO MIN HECKONBKH AKKYMYJIATOPOB, AKKYMy- 
AHPYIOlNNS SHAyeHUA PyHEUMÜ MIA HEKOTOPBIX SHAYCHNÄ APryMeHTa U HHTEPHONALMOHHBIX 
hakrtopoB. Ms aTux sHayeunt PyHRumä O060pasyWmTcH IyTeM IPOCTOü MIN CIIOIKHOÄ Mexa- 
HUYeCKOU UHTEPIIOAHALUU BHAYEeHNA PYHKUHÜ IK JIOÖBIX APTyMeHToB. Kparko YKasblBalmTcH 
TeXHHYeCKUe BO3MOKHOCTU BEIYUCHEeHUH MH MECEeNyeTcH TOYHOCTB, KOTOPAA MOIKeT ÖBITb 
nocruruyra. B cakrımycHne 3TOT UHTEePIONAUMOHHBIH CIIOCOO CPABHHBAETCH CAPYTHMH MeETO- 
MaMmn. 


In dem Vortragsreferat ‚Die Funktionsrechenmaschine‘‘t) wurde angedeutet, daß sich 
der Rechenbereich einer numerischen Rechenmaschine über die vier Grundrechnungsarten 
hinaus dadurch auf Funktionen ausdehnen läßt, daß die Maschine mit einem Speicherwerk 
versehen wird, in welchem für grob abgestufte Argumente genäherte Funktionswerte und Inter- 
polationsfaktoren gespeichert sind, aus denen die Funktionswerte für beliebige Argumente 
durch maschinelle Interpolation gebildet werden. Die Interpolation kann einstufig (linear) 
oder mehrstufig erfolgen. Die Argumentabstände können gleich oder verschieden sein. Die 
sich daraus ergebenden technischen Lösungsmöglichkeiten und Interpolationsfehler sollen nach- 
stehend näher beschrieben werden. 


I. Einstufige Interpolation 
I. 1. Prinzip der einstufigen Interpolation 
Um die Funktion f(x) im Bereich ©, <x< x,, darzustellen, unterteilen wir den Bereich 


in m Intervalle, die zwischen den ‚‚Grundargumenten“ &,, &1, &g -: : Lyss int is: : Im—1, Im 
liegen. Innerhalb jedes einzelnen Intervalls ©,< = < &, + ı ersetzen wir f(x) durch eine Gerade 
=. tb, Zn) = mtb Az....:.... ll) 


Dabei werden die Faktoren «a, und 5, jeweils so gewählt, daß g,(x) in dem betreffenden Inter- 
vall möglichst wenig von f(x) abweicht. Speichern wir die Werte a,, b, für jedes einzelne Inter- 
vallin einer Rechenmaschine und versehen die Maschine mit einer Einrichtung zum Übertragen 
der Werte aus dem Speicherwerk in die Rechenmaschine, so läßt sich der Näherungswert g,(2) 
für jedes beliebige Argument &—= x,-+4Ax innerhalb des gespeicherten Bereichs dadurch 
berechnen, daß das Speicherwerk entsprechend dem Grundargument x, eingestellt wird und 
daß beispielsweise der Wert a, in das Resultatwerk, der Wert b, in das Einstellwerk der Rechen- 


j _ maschine übertragen wird. Multiplizieren wir anschließend 5, mit dem „Restargument“ Az, 


25 


so steht im Resultatwerk der gesuchte Näherungswert von f(@). 


2) Z. angew. Math. Mech. 30 (1950), S. 294—295. 
2 21 
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I. 2. Technische Lösung der einstufigen Interpolation 
Die einstufige Interpolation läßt sich technisch einfach lösen. Nach der Patentanmeldung 
von 1941 kann z.B. die Funktion sin x wiein Bild 1 dargestellt gebildet werden. In dem Speicher- 
werk 9 sind Näherungswerte von sin & (Grundwerte a,) für die Grundargumente 0°, 1°, 2°... 89° 
und die zugehörigen Tafeldifferenzen für 1 Minute 5, ...-Ög; in 90 Blechen 8 durch verschieden 
tiefe Einschnitte dargestellt. Es soll nun z.B. 


sin 35T 247 Frag; Das > DAT WR er 


berechnet werden. Zu diesem Zweck wird das Speicherwerk 9 durch Drehen, von Hand oder 
automatisch, auf das Grundargument 35° eingestellt, so daß die zu 35° gehörigen Bleche 10 
und Il in die Abgreiflage kommen. Hierauf werden die Zahnstangen 12 (z. B. durch Federkraft) 
aus der Nullstellung bis zum Anschlag an Blech 10 verschoben, wodurch der Wert a, über die 
Zahnräder 13 in das Resultatwerk 14 übertragen wird. Gleichzeitig werden die Zahnstangen 15 
aus der Nullstellung bis zum Anschlag an Blech 11 bewegt. Dadurch wird die Tafeldifferenz b;; 
über die Zahnräder 16 und die Einstellsegmente 17 in die Schaltklinkenräder des Einstellwerks 18 
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einer normalen Schaltklinkenmaschine übertragen. Nunmehr wird das Resultatwerk 14 unter 
das Einstellwerk 18 geschoben und die in 18 eingestellte Tafeldifferenz mit der Anzahl der Minu- 
ten, also mit 24 in der üblichen Weise multipliziert. Im Resultatwerk 74 steht dann der gesuchte 
Funktionswert. 

Die Konstruktion läßt sich sehr vielseitig variieren, so daß jede Maschinentype mit einer 
oder mehreren ‚„‚mechanischen Funktionstafeln‘, die zudem noch auswechselbar sind, ausgerüstet 
werden kann. Die Grundwerte und Tafeldifferenzen können auch durch Lochkarten, elektrische 
Drehwähler, Relais usw. gespeichert werden. Die Übertragung aus dem Speicherwerk in die 
Rechenmaschine kann mechanisch, elektrisch-mechanisch oder rein elektrisch erfolgen. Die 
Funktionswerte können voll- oder halbautomatisch oder durch Handbedienung gebildet werden. 
Verwenden wir eine Doppelrechenmaschine, so können wir in einem Rechengang zwei ver- 
schiedene Funktionen, z. B. sin x und cos &, berechnen. Auch lassen sich etwaige Vorzeichen- 
und Quadrantenregeln leicht automatisieren. 


Die Übertragungseinrichtungen können auch mit dem Umdrehungszählwerk U und dem 


Resultatwerk R bzw. mit U und dem Einstellwerk E zusammenwirken. Im ersten Fall wird 

b„in U, a, in Rund Axin E eingestellt, worauf nach einem Multiplikationsgang g,(x) in R 

erscheint. Im zweiten Fall wird —a, in U, 1/b, in E und Axzin R eingestellt. Dann erscheint 

nach einem Divisionsgang —g,(x) in U. Verwenden wir nur eine einzige Übertragungseinrich- 

lung, die mit E, R oder U zusammenwirkt, so entsteht g,(x) auf folgende Weise: 

1. Übertragung nach E: Übertragen von a, nach E und von dort durch Multiplikation mit 
Eins nach R. Übertragen von b, nach EZ. Multiplikation mit. Ax. In R steht g,(«). 


2. Übertragung nach R: Übertragen von b, nach Rund von dort durch Rückübertragung nach Z. 
Übertragen von a, nach R. Einstellen von Axin U. Multiplikationsgang, in R steht (8). 


3. Übertragung nach U: Übertragen von a, nach U und von dort durch Multiplikation mil 
Eins nach R. Übertragen von b, nach U. Einstellen von Axin E. Multiplikationsgang, in 
R steht g,(2). 


Werden die gespeicherten Grundwerte a, ziffernmäßig angezeigt, so kann die mechanische 
Funktionstafel auch dazu benutzt werden, zu einem gegebenen Funktionswert f(@,+4Ax) das 
zugehörige Argument 2,—+Ax zu ermitteln. Dies kann z.B. folgendermaßen geschehen: 
f(z.-+-4Ax) stehe nach einer vorausgegangenen Rechnung in R oder wird dort eingestellt. Nun 
wird das Speicherwerk so eingestellt, daß der nächst kleinere oder größere Wert a, in die Ab- 
greiflage kommt. Dann wird a, nach E übertragen und von f(&, + Ax)in Rabgezogen. Anschlie- 
ßend wird b, nach E übertragen. Dividieren wir jetzt die Differenz f(&, + Ax) — a, in R durch 
b,, so erscheint in U Az, das zu dem am Speicherwerk abgelesenen Grundargument z, addiert 
wird. 

Wir können auch nur die Grundwerte speichern und die Tafeldifferenzen mit der Maschine 
berechnen, z.B. auf folgende Art: Wir übertragen a, nach E und von dort durch einen Addi- 
tionsgang nach R. Dann wird a, nach E übertragen und von «,+1in Rabgezogen. Die Differenz 
@n +1 —q, wird mit Hilfe einer Rückübertragungseinrichtung aus R nach E übertragen und an- 
schließend mit Az/h 2) multipliziert. Schließlich wird a, nochmals nach E übertragen und von dort 
durch einen Additionsgang nach R gebracht. Dieses Verfahren hat den Vorteil, daß das Volumen 


der mechanischen Funktionstafel nur halb so groß ist wie bei Speicherung von Grundwert + Ta-. 


feldifferenz. Nachteilig ist jedoch, daß sich Ax/h von Ax nur durch die Kommastellung unter- 
scheiden darf, daß wir 3 Übertragungen, 2 Additionen, 1 Subtraktion, 1 Multiplikation und 
1 Rückübertragung brauchen gegenüber 2 Übertragungen, 1 Addition, 1 Multiplikation und 
0 Rückübertragungen, wenn die Tafeldifferenzen mit gespeichert werden. Die Rechnung geht 
deshalb im letzteren Falle schneller und einfacher und ist auch anwendbar, wenn Az/h ziffern- 
mäßig nicht mit Ax übereinstimmt. 
I. 3. Genauigkeit der einstufigen Interpolation 

Bei der üblichen linearen Interpolation wird f(x) im Intervall «, <a,+4dx2 <a, h 

durch die Sehne ersetzt. Dadurch entsteht ein Interpolationsfehler, der seinen Maximalbetrag 


em; 3 
LE 

etwa in der Intervallmitte erreicht. Dabei ist 2, <&£< x,-+ h. Ersetzen wir die Sehne durch 
eine vermittelnde Gerade, so geht der maximale Interpolationsfehler auf etwa die Hälfte zu- 
rück. Er tritt dann am Anfang, am Ende und etwa in der Mitte des Intervalls auf und beträgt 


lünael = 178) - RR N 


2) k — Intervallbreite, 


« 
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Der Interpolationsfehler ist also abhängig von der Intervallbreite A und der zweiten Ableitung 
der gespeicherten Funktion. 

Wenden wir Gl. (4) auf die wichtigsten trigonometrischen Funktionen an und setzen & =, 
so erhalten wir für neue Teilung mit je 100 Grundwerten und 100 Tafeldifferenzen im Speicher- 
werk bei gleichen Argumentabständen die in Tafel 1 zusammengestellten Interpolationsfehler. 


Tafel 1 


Zusammenstellung der maximalen Interpolationsfehler der wichtigsten trigonometrischen 
Funktionen bei einstufiger Interpolation und konstanten Intervallbreiten 


Bereich h? Max | Umazx | 
' h le. |f" 
si & "Im 4 Omas | IR Betrag | Stelle 
sin © 08 + +» 1008 18 | 1,54 - 105 - sin x 1,54 105 = 1008 
cos 08 » - » 1008 18 1,54 - 105 - cos 1,54 - 10 = (08 
tg & 08 +. - 508 0,58 0,77. 10-5. BE 1,54 + 10-5 x — 508 
cos?x 
etg & 508 + + + 1008 0,58 0,77 . 10-5. 18 ® 1,54 - 10-5 x — 508 
sin? x 
arc sin « 0--+1 0,01 3,9900 2. E1 — 1,38 3) are sin &— 1008 
(08 - - » 1008) (1-29? 
arc cos © 1’--0 0,01 3,98eC » = 2 1,38 3) arc cos 2 —= 08 
(08 +» - 1008) (1-22) 
x 
arctg x 0-7 0,01 7,96°° . —_— — 2,58°° arctg x = 348 
2 (08 + + - 508) der | 
are ctg © 1---0 0,01 7.96c8 . T 2,58ee are ctg = = 668 
| (508---1008) ° a EN 


Bei tg x empfiehlt es sich, nur den Bereich 08... 508 darzustellen, weil f"(x) für z— 1008 
gegen Unendlich geht. Für 508 <z<1008 berechnen wir tg (100°—x) und erhalten dann tgzaus 
tgx = 1/tg (1008—x). Entsprechendes gilt für ctg«. Bei arcsin x und arc cos x wirken sich 
gleiche Argumentabstände in der Nähe von = 1 sehr ungünstig aus. Hier kann jedoch, wie 
später noch gezeigt wird, durch eine geeignete Verteilung der Grundargumente die Genauigkeit 
wesentlich gesteigert werden. Dagegen sind die Funktionen arc tgx und arc ctg& für die ein- 
stufige Interpolation sehr geeignet. Der maximale Interpolationsfehler beträgt hier nur 2,58°, 
Dabei läßt sich durch einfache konstruktive Maßnahmen erreichen, daß mit einer mechanischen 
arc tg- bzw. arc ctg-Tafel mit einem unmittelbaren Rechenbereich von nur 50% jeder beliebige 
Winkel berechnet werden kann, ohne daß der Rechner Vorzeichen- oder Quadrantenregeln be- 
achten muß. 

Für Funktionstafeln in alter Teilung erhalten wir für sinz und cos& für = 1° und für 
tgx und ctgz für h = 0,5° Interpolationsfehler, die 1,23 mal größer sind als die entsprechenden 
Fehler bei neuer Teilung. Für arcsin® und arc cosz wird Max|v„..| = 1,2°, für arctgx und 

„arc ctgz wird Max |v,.a.| = 0,84". 

In Tafel 1 sind die Intervallschritte innerhalb der einzelnen Funktionen gleich groß und die 
Grundargumente runde Werte, wodurch sich das Restargument Az stets auf dieselbe Art und 
Weise ergibt, was für einen vollautomatischen Rechengang von Vorteil ist. Wird das Speicher- 
werk von Hand eingestellt, so können die Intervallschritte auch ungleich sein. Wählen wir A 
stets so, daß in jedem Intervall der zulässige Interpolationsfehler gerade erreicht aber nicht über- 
schritten wird, dann erhalten wir mit einem Minimum an gespeicherten Werten die gewünschte 
Genauigkeit. Damit sich jedoch Ax noch leicht im Kopf rechnen läßt, werden zweckmäßig noch 
einige. Intervalle eingeschoben, so daß die Grundargumente nicht gar zu unrunde Werte anneh- 
men. Stufen wir z.B. die Grundargumente für sin x wie in Tafel 2 angegeben ab, so erhalten wir 


über den ganzen Bereich mit nur 134 Grundwerten und 134 Tafeldifferenzen eine gleichbleibende 


Genauigkeit von etwa 5 - 10”, d. h. dieselbe Genauigkeit wie bei einer gedruckten fünfstelligen 
Sinustafel. 
Tafel 2 


ZweckmäßigeIntervalleinteilung einer fünfstelligen mechanischen Sinustafel mit ungleichen 
Intervallbreiten und einstufiger Interpolation 


Bereich 2 7... 0,0—4,0—10,0—11,4—12,7—15,1— 19,5 — 24,5 —29,9—40,3—55,0 — 1008 
Intervallbreite . . ...r 2,0. 1,5. 14 1,3-771,227,1 7,0720 DEE BEE RE 


Zahl der Intervalle... 2 all eı u 1. CD SEAeEEe BE 
®) Unmittelbar als Differenz zwischen Kurve und Interpolationsgerade berechnet. =, 


Ne 
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Wird die mechanische Funktionstafe] mit einstufiger Interpolation dazu benutzt, um aus 
einem gegebenen Funktionswert das zugehörige Argument zu bestimmen, so folgt aus 
RS sm 
mt A)=n tb da tee EN HR TELON 
WOLLT 27 ist, 
_ Fan + Ar) — a, m f'(&) 
RT ale ER et (6). 
Der Fehler von Az und damit von & ist also 
2 " 2 "„ e 
ae MO, at M| 


- =&E'— 1, ——— 0.2... 0... un 
| 16° 78) TE nA] m 
wo 0<d<1l und O<n<{1 ist. Der Maximalfehler ist somit 


_R |\SÜ+Bh| PR fa) 
deln, Na Eur ENT ET LIEST 
Für f(x) =sinx wird 
[Alma = tg N RE RN RAR) 


Wenden wir Gl. (9) auf eine mechanische Sinustafel entsprechend Tafel 2 an, so erhalten wir für 
- ) den Bereich 55 < x< 1008 |d&| ar 7 3,5°° tg @,. Das entspricht für den angegebenen Bereich 
etwa dem Maximalfehler einer gedruckten fünfstelligen Sinustafel |d&|nas 7 3,2°° » sec ®. 


II. Mehrstufige Interpolation 
1I. 1. Prinzip der mehrstufigen Interpolation 


Die Genauigkeit der einstufigen Interpolation ist begrenzt durch die Anzahl der Werte, die 
im Speicherwerk untergebracht werden können, letzten Endes also vom zulässigen Volumen des 


Speicherwerkes. Wollen wir die Genauigkeit erheblich steigern, ohne das Speichervolumen 


wesentlich zu vergrößern, so müssen wir mehrstufig interpolieren. Zu diesem Zweck unterteilen 
wir den Bereich 2,<z< x, in welchem die Funktion f(x) maschinell berechnet werden soll, 
in m Intervalle gleicher oder verschiedener Breite h. Innerhalb eines jeden einzelnen Intervalls 
%,<&<®,+ h ersetzen wir f(z) durch eine Näherungsfunktion 


Den rd m AR + AP: RR "ST He (10). 
Dabei ist wieder 

2 Se a ee SR a Re (11). 

Schreiben wir Gl. (10) ir der Form 
In (an + Ax) = + 4 (&,+ A% (0, +42: ))) RN A (12), 
so kann sie mit einer Rechenmaschine, in welcher die Faktoren «a,, b,, 0, d„ . . . gespeichert sind, 
ohne nennenswerte technische Schwierigkeiten halb- oder vollautomatisch gelöst werden, wie 

L) nachstehend gezeigt wird. 


II. 2. Technische Lösung der mehrstufigen Interpolation 


In Bild 2 ist ein Ausführungsbeispiel einer Funktionsrechenmaschine mit mehrstufiger 
Interpolation schematisch dargestellt. Eine Rechenmaschine mit dem Einstellwerk E, dem 
Resultatwerk R, dem Umdrehungszählwerk U und 
einer Rückübertragungseinrichtung von R nach E ist 
zusätzlich mit einem Einstellwerk A für x,, einem 
Speicherwerk S, für Az, einem Speicherwerk 8, für 
die Faktoren a,, b,, 6,, d, . ... und Übertragungseinrich- 

tungen von S,nach U und von S, nach Eversehen. Damit 
u läßt sich Gl. (12) folgendermaßen lösen. Wir stellen z, in 
A und Azin S,ein. Durch das Einstellen von «, wird 


das Speicherwerk 8, so verstellt, daß die Faktoren a,, b,, 

0 Ip. . . indie Abgreiflage kommen. Anschließend wird Mi ER 
dd, nach E und As nach U Eu Rasen. Nasa äihem E eohenmaschine nıft mehretufiger Interpolation, 
» Multiplikationsgang steht in x *d,. Dann wir 
3P ‘c„ nach E übertragen und durch Addition (positiv oder negativ) nach R gebracht. Die Summe 
+ 4&:d, wird nach E rückübertragen und mit As multipliziert. Nunmehr wird b, nach E 


j ER 


nd von dort durch einen Additionsgang nach R gebracht. Die Summe db, + 4x (+42 (d„+*)) 


ar 


."Multiplizieren wir die Klammern aus, ordnen nach Potenzen von Ax und ersetzen k durch h/p, 
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wird wieder von R nach E rückübertragen und mit Ax multipliziert. Schließlich wird noch «a, 
nach E und von dort durch einen Additionsgang nach R gebracht, so daß dort der gesuchte Wert 
9,(@, + Aa) zf(@,+ Ax) steht. Es wiederholt sich also immer wieder dasselbe Rechenspiel: 
Multiplikation mit Ax — Übertragung aus 8, — Addition — Rückübertragung von Rnach E — 
Multiplikation mit Ax usw. Auf die Zablröichen technischen Variationsmöglichkeiten kann hier 
nicht eingegangen werden. Erwähnt sei nur, daß bei zweistufiger Interpolation mit einer Hand- 
rechenmaschine auf ein besonderes Speicherwerk für Ax verzichtet werden kann, indem bei der 
ersten Multiplikation Ax positiv oder negativ in das Umdrehungszählwerk eingekurbelt und bei 
der zweilen Multiplikation wieder herausgekurbelt wird. 


Il. 3. Genauigkeit der mehrstufigen Interpolation 
Die Fakloren a,, b,, 6, . .. sind in jedem Intervall so zu bestimmen, daß g, (x) von f(x) mög- 
lichst wenig abweicht. Dies können wir z.B. folgendermaßen erreichen: Wir zerlegen jedes 
Intervall ©, << x,-+- hin p gleiche Teilintervalle von der Breite h/p= k und fordern, daß 
9„(x) an den Stellen %,, hl k,%,-+2°k...2%,-4p:% genau mit den Funktionswerten 
fa), fan +1:%, fu + 2°k)... (+ D:%) übereinstimmt. Diese Forderung wird erfüllt 
von der Newtonschen Interpolationsformel 


9n(%) —=f( %,) Ze (7 Ln a + (2 — %,) (x ro 13 k) y a 


| (13). 
+1 de). 


Die Differenzen Af(z,) usw. ergeben sich aus dem Differenzenschema 


urn ®.) 
| Af(@) 
„+1:k a+1:.k) | Pf.) 
Af(a +1 k) Af@,) 
%+2:% f.+2:% An 1 Mr 
Ale 2: AN a N ee ne 


%+3:% BR NEN 2 N a Ar Re 


%,, = D. k fan H Pu k) 
Setzen wir ©— x, = Az, so erhalten wir aus (13) ; 
I) =f@) + de As-(As—1 Be | a 


+ + 4a ch). A a 


so wird 


A Ar ee Arsen) u ER 


14a. ee Ale). nn N | | - » a 
ai \ 


Gl. (14) gilt unmittelbar für vierstufige Interpolation, wenn wir p = 4 Selen Für ee In f 
terpolation sind die Glieder mit Differenzen der Ordnung p+ 1 und höher“ egzulassen. Das 2 “ [A 
Glied f : die Koeffizienten von Az, Aa? .. .in Gl. (14) ERIPRe che den re As Da» 6 re 

von Gl. (1 . 


Der Fehler, den wir bei Anvenäiiag der Newtanen, erplation rn 


„fa, - CE ZmWe-m-ih... 


ER ETF 
Fa Be Br 


Bee VA Der 1 a a neh hal End En Be a) =. Zar Kuh ES 
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wo £&,eine Stelle zwischen &, und«,+p: k= 07 —+ h bedeutet. Setzen wir k= h/pund &— 2, 
=Ar=u:h wo 0<su<s[l ist, so erhalten wir 
WE u(u — 1/p) (u — 2/p)... 


5 o 7 ni (u Es pp) pri :fw +1) (&,) — U, hp+i fo +1) (£,) (15). 


In Bild 3 ist U, für ein- bis vierstufige Inter- a er Er ee. = 
polation dargestellt. Wir entnehmen daraus \ I u 
[U ma21 1,25 ..10°2.an der ‚Stelle u = 0,5, | i 


|U, nar| = 8,02 : 103 an der Stelle u 0,211 
und «= 0,789, 
[Ugmaz| = 5,14 : 107% an der Stelle u = 0,127 
und u 0,873, 
| Uymar\ = 2,96 : 10° an der Stelle «= 0,089 
und % = 0,911. 
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Die maximalen Fehler. einer p-stufigen Inter- 
polation sind somit: 


d 1 [Omas] —(1,25 EIUR ae ar (&)); 
wi ne 18:02 1078 73>f77 (£&,)|, (16). 
P=3 |Osmas| 514 - 10. HM. fIV&)], 
p — 4 |Umaz| —=1|2,96 e 105 - RP - RN. (&)] 
Sie lassen sich durch ‚‚vermittelnde Interpola- 


tion‘ verkleinern. Zu diesem Zweck schreiben 
wir Gl. (15) in der Form 


N den a 5) 2” 
en U,) mr ferne) di] ( " 


und bestimmen U, so, daß der zugehörige Teil- 

fehler v, durch Ändern der Koeffizienten von 

9n(#) berücksichtigt werden kann, und das 

Maximum des verbleibenden Restfehlers 20 02 04 06 8 uU m 
|%p a 5 |% a 

ist. In Tafel 3 sind die erreichbaren %p maz und 


zugehörigen Verbesserungen v, von g,(z) zusammengestellt. Die Funktionen U, sind in Bild 3 
dargestellt. Dort ist auch ersichtlich, wo |U,—UD,| und damit |v,| ein Maximum "wird, 


Bild3. Graphische Darstellung von U,, und Do 


Tafel 3 
Zusammenstellung der maximalen Fehler |%, man! bei vermittelnder Interpolation und der- 
Verbesserungen v, von gn (%) 


| % Er Hrönaalı 
1 | —6,25 10-2. R2- f"(£,) 6» 1072. R2- If" (£))] 
2 | 5-10#. nf (&,) — Aw 1072. 027 (&,) 5.109. 9.1’) 
3 | 3-10. A-fIV(&,) + Au 12-104» RafIV(&,)— Aa? 12-104 HB fIvVlE)| 3-10-4- m |fIVE,)) 
4 | 1,5-105- 15. fv(&) — 4a 3-10 1 fV(&,) 1,5 + 10-5. 75- |fV + (&,)) 
Setzen wir Rz Ber. a U N a ER REG :)" 
so wird " l 
Na A ET Fe PR RE RE) 


In Bild 4 ist U, in Abhängigkeit von h und p dargestellt. Mit Hilfe dieses Diagramms und der 


| (p-+ 1)-ten Ableitung der darzustellenden Funktion f(x) können wir für einen vorgegebenen Ab- 


stand A der Grundargumente die maximalen Interpolationsfehler und für eine vorgeschriebene 
j ‚Genauigkeit das erforderliche h ermitteln. So erhalten wir z.B. für f(x) = sin x und h = 1! 


[öfmas| = 1,5103 sin &ı|»  1V%maz| = 11,8 104 - sin &l, | 0 
[95 maz | = [1,9 : 10°® - cos &l» [0" max | = 11,410 cos &|. I (20). 


N 
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Durch Speichern von 2 - 100, 3 : 100, 4 - 100 bzw. 5 : 100 Interpolationsfaktoren kann also etwa 
die Genauigkeit einer 5-, 8-, 11- bzw. 14-stelligen Sinustafel erreicht werden. Soll der Interpola- 


92 10-2 tionsfehler |vy:.:| nicht überschreiten, 
so erhalten wir die zulässige Intervall- 

“+ 10% breite A, indem wir in das Diagramm, 
Bild 4, mit p und 

Up De = [op 2uil \f® ir en Mean 

0° ws eingehen. So ergibt sich z.B. für 


Fa) sin z, \f? RN (Elan > 1, 
op zur == 10, D% Mhz 107 


7072 70° 72 
für zweistufige Interpolation 
hy NT 0,88 
0 417076 und für dreistufige Interpolation 
— . h; BUDAT 58, 
20 FEAR RT. 20 Wir können also die Zahl der zu spei- 
70° 70° n x : dar 
000° 0002 000% or So a Br chernden Werte ohne Grenauigkeits- 
verlust erheblich reduzieren, wenn wir 
079 029 EL 29 39 59 die Zahl der Interpolationsstufen ver- 
mehren. Innerhalb jeder Interpola- 
107 320 SORR EN FE an DES tionsstufe kann die Zahl der zu spei- 
Bild4. Graphische Darstellung des Faktors U, zur Berechnung des maxi- chernden Werte noch dadurch verklei- 
malen Interpolationsfehlers nach der Gleichung op maz\=| U en En]. nert werden, daß h nicht konstant son- 


dern jeweils so gewählt wird, daß der 
zulässige Interpolationsfehler etwa erreicht aber nicht überschritten wird. Dies ist der Fall, 


D+ ll ——— 
wenn h proportional a (£,) gemacht wird. 


Die ein- und zweistufige Interpolation kommt in Frage für Rechenmaschinen mit kleiner 
Rechengeschwindigkeit. Für schnellrechnende Maschinen, bei denen die Zahl der Rechenope- 
rationen nicht sehr ins Gewicht fällt, ist dagegen eine 3- oder mehrstufige Interpolation vorteil- 
hafter, weil dadurch das Volumen des Speicherwerks bei gleichbleibender Genauigkeit erheblich 
verkleinert werden kann. 


Ill. Vergleich des Interpolationsverfahrens mit anderen Verfahren 
Funktionswerte lassen sich auch ohne Interpolation maschinell berechnen. So können wir 
z.B. bei den trigonometrischen Funktionen auf die Additionstheoreme zurückgreifen. Speichern 
wir z.B. die Werte sin », und cos &, von 18 zu 18 und cos Axund sin Azim Bereich 0 < Ax < 100° 
von 1° zu 1°, so können wir für auf volle Minuten abgerundete Zwischenargumente sin & berech- 
nen aus 
2 sin (2, + Ax)= sing, cos A& -{- cos 2, :sin Aua. 2.02 2.2.0212). 


Dieses Verfahren versagt jedoch, sobald die Funktion für beliebige Ax mit großer Genauigkeit 
berechnet werden soll. 

Denselben Nachteil hat der Vorschlag, jeden einzelnen Funktionswert unmittelbar zu ver- 
körpern, also z.B. von der Funktion sin x die Werte für x = 0°, 1°, 2°... 100 00° zu speichern. 
Ein weiterer Nachteil dieses Verfahrens ist, daß es schwierig ist, 10000 Werte auf kleinem Raum 
zu speichern, und eine derartig umfangreiche mechanische Funktionstafel schnell ‚aufzu- 
schlagen‘. \ 


Das andere Extrem ist, die Funktion (x), falls dies möglich ist, in eine Reihe 
S)=a+b-a te R+d- RP 4... ...02.2.22.(2) 


zu entwickeln, die für den ganzen Bereich, in dem f(z) dargestellt werden soll, gültig ist. Schrei- 
ben wir Gl. (22) in der Form 


f)=a+z(lb+sa(c+z(d Er N. 2 anno): 


so haben wir denselben Aufbau wie bei den Näherungsfunktionen g,(x) für mehrstufige Inter- 


polation Gl. (12). Wir können deshalb Gl. (23) entsprechend II. 2. maschinell lösen, wenn wir I 
die Werte a,b, c... speichern. Wir kommen in diesem Falle mit einem Minimum an gespeicher- 
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ten Werten aus, müssen dafür aber sehr viele Rechenoperationen ausführen um die gewünschte 
Genauigkeit zu erreichen. Für sin x erhalten wir z.B. wenn z in Neugrad gegeben ist, 


nee} 


Wollen wir nach (24) sin x auf 7 Stellen genau erhalten, so müssen wir für x = 1008 bis zum Glied 
mit ol®rechnen. Hierzu sind 13 Multiplikationen, 12 Rückübertragungen und 6 Additionen bzw. 
Subtraktionen notwendig. Dieses Verfahren eignet sich also nur für sehr schnell rechnende 
(z.B. elektronische) Rechenmaschinen. 

Das Interpolationsverfahren nimmt also eine Mittelstellung ein zwischen der unmittelbaren 
Verkörperung jedes einzelnen Funktionswertes, bei der am meisten Werte zu speichern sind und 
am wenigsten zu rechnen ist und dem Reihenentwicklungsverfahren, bei dem am wenigsten 
Werte zu speichern sind, dafür aber am meisten zu rechnen ist. 


Berechnung der turbulenten Strömung längs der ebenen Platte 
Von W. Szablewskı in Berlin 


Die Differentialgleichungen der turbulenten Grenzschicht an der längsgeströmten ebenen Platte werden 
unter Zugrundelegung des Prandtlschen Schubspannungsansalzes von der Wand her in zwei Schritten 
integriert. Die so gewonnene Lösung dient für den übrigen Bereich der Grenzschicht als Approwimation. 
Bei passender Wahl der empirischen. Koeffizienten der Theorie erhält man gute Übereinstimmung mit den 
experimentellen Ergebnissen. 

For a viscous fluid flowing along a plane plate, the differential equations of the turbulent boundary 
layer are integrated, by use of Prandils relation for Ihe shearing stresses, in iwo steps starting from the wall. 
The solution obtained is taken as an approximation for Ihe remaining, part of the boundary layer. Good 
conformity with Ihe experimental dates is reached by appropriate values of the empirical coefficients oceuring 
in the theory. 

Pour la couche limite d’un fluide coulant le long d’une surface plane les Equations diff£rentielles son 
inlögrees en deux pas a parlir dela paroi; on y a fait usage de la relation de Prandil pour l’effort de cisaille 
ment. La solution tellement obtenue peut servir comme approwimation pour le reste de la couche limite, 
En choisant convenablement les coefficients empiriques de la theorie on obtient bonne correspondance avec 
les resultats ewperimentauw. 

udeperumansHusie ypaBHeHuA BUXPeBOTO TOTPAHHYHOTO CIAIOA TIOTOKA BAONB TNIOCKOCTH 
HHTerpnpymTcaA pn moMomm dopmyısı Ilpanıraa yıa kacaTebHbIX HANPpsIKeHnÜ B Ba ara 
B HAIIPABJICHHH, HEPHEHAHKYIIAPHOM K OÖTekaemoh IIocKOcTu. llosyyeHHoe pemenHne cHy3KuT 
TIpnönmkeHHoä PopMyAloH MIA OcTamımefich yacTu MOTPaHnyHOoTO cı1oA. Ilpıu monxonsmem 
nonÖope 3MIIHPHYeCKUX KOSPHLMEHTOB TEOPHH MNONYYAeTCcH XOopomee COBHANEHHE C MAHHEIMH 


OLEITOB. 
Bezeiehnungen 
o Dichte. UN Tr i : Pie: 
Kanne y= — Dimensionslose Geschwindigkeit. 
u Zähigkeit. %% 


RR RRKUE h 
WIRD lalsche, Zalugbeit,. x *7 Dimensionsloser Wandabstand. 
v 


u, v Geschwindigkeitskomponenten. 


u, Anströmgeschwindigkeit. " Dimensionslose Dicke der laminaren 
v, Querkomponente der Strömung am u. 0 | Unterschicht. 
Rande der Grenzschicht. v 


M) =, Dimensionslose Geschwindigkeit. I Mischungsweg. 


x, u Koeffizienten der Mischungsweg- 


x Abstand von Plattenvorderkante. \ 1 y N 
y Wandabstand. verteilung „— x, Tu (4 
6 Grenzschichtdicke. dir 
Ru — 
Yy . ö 


1-5 Relativer Wandabstand. 


Reynoldssche Zahlen. 


tr Schubspannung. 
to Wandschubspannung. 
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2a. 1 i Koeffizient der Ge- 
Tu Fi I+ 2) schwindigkeitsver- 
= f Örtliche Widerstandszahl er 
= rtliche Widerstandszahl. 
/ eo r w— Lei {Inn —Cn} + const. 
2 
5 | =] 1—o)dn Verdrängungsdicke der Strö- 
W= [ x,dx Plattenwiderstand. Bu: 
0 
“ Ä eo ee 
0, Gesamtwiderstandszahl. 92 RR 2 3 
2 | ) 
» U, % 4A — 4 BR a,el +C—x#%o . 
107 C ’ Ableitung nach «. 
DE 27 > Wenig verschieden von. 


I, Einleitung 
Der vorliegenden Untersuchung der turbulenten Strömung längs der ebenen Platte wird 
für die scheinbare, durch den Vorgang der turbulenten Mischung erzeugte, Schubspannung — wie 


in [1] — der PrandtlIsche Schubspannungsansatz 
Pe ou ou 
5 jaylay 
mit dem Mischungsweg / zugrundegelegt. Für den Mischungsweg machen wir den Ansatz 
I Re & Habe A A 
BEN HR BE G sc } 5 
"5 Fu 5 -|- (ö Grenzschichtdicke) 


Die Untersuchung der Strömungsvorgänge an der ebenen Platte erweist sich erheblich schwieriger 
als die Untersuchung der .Rohrströmung [1]. Zunächst einmal ist die Reynoldssche Zahl 


Re,= - (%, Geschwindigkeit der Anströmung, x Abstand von der Plattenvorderkante) längs 


der Platte mit wachsender Entfernung von der Vorderkante variabel, so daß wir also mit von 
Strömungsquerschnitt zu Strömungsquerschnitt variierenden Geschwindigkeitsprofilen zu 
rechnen haben. Entsprechend hätten wir im Mischungswegansatz die Koeffizienten infolge ihrer 
Abhängigkeit von Zähigkeitseinflüssen bzw. der Reynoldsschen Zahl als längs der Platte 
variable Größen aufzufassen. Mit Bezug darauf aber, daß für große Reynoldssche Zahlen der 
Zähigkeitseinfluß ein sehr geringer wird (s. die entsprechenden Darlegungen in [1]), werden wir 
bei Beschränkung auf große Re-Zahlen die Koeffizienten x und „ als Konstanten ansehen. Das- 


selbe gilt für die in den Rechnungen auftretende dimensionslose Dicke en der laminaren 


1/To : 
Unterschicht | — y= . Des weiteren kommt als die mathematische Untersuchung sehr er- 


schwerend das Auftreten der (bei der Rohrströmung fehlenden) Trägheitsglieder hinzu. 


Wir werden eine stufenweise Integration der Bewegungsdifferentialgleichung von der Wand 
her — wie bei der turbulenten Rohrströmung — vollziehen. Im ersten Schritt werden wir die 
Integration über die Jaminare Unterschicht, den sich anschließenden Bereich gleicher Größen- 
ordnung von laminarer und turbulenter Reibung bis in den Bereich der (überwiegend) turbulenten 
Reibung hinein vollziehen. Bei diesem ersten Schritt ist der Integrationsprozeß im Bereich der 
turbulenten Reibung beschränkt auf ein Gebiet, für das 1 »-xy (=> wenig verschieden von). Die 


Trägheitsglieder sind — wie die Untersuchung zeigen wird — im Bereich des ersten Integrations- 


schrittes von kleinerer Größenordnung und brauchen nicht berücksichtigt zu werden. 

In einem zweiten Schritt werden wir dann die erhaltene Lösung fortsetzen. Da wir bereits 
beim ersten Integrationsschritt bis in den Bereich der (überwiegend) turbulenten Reibung gelangt 
sind, können wir beim zweiten Schritt nunmehr die laminare Reibung vernachlässigen. Im 
Bereich der turbulenten Reibung dehnen wir den Integrationsprozeß auf ein Gebiet aus, das 


Ä ‘ I A zu 
gekennzeichnet ist durch < =* y + „(2) . Insofern wir die turbulente Reibung auch in 


Ö () ) 
zweiter Größenordnung berücksichtigen, wird es sich als notwendig erweisen, nunmehr auch die 
Trägheitsglieder zu berücksichtigen, die dieselbe Größenordnung besitzen. 


Die so gewonnene Funktion dient im übrigen Bereich der Grenzschicht als Approximation. = 2 


Grenzschichtdicke ö und Wandschubspannung z,, die in der Funktion als Parameter eye 
werden dann mittels einer Randbedingung und des Impulssatzes berechnet. EREN - 


a k 


3 
“ 
3 


{ i Bu RR ; e. | = ‚u ER 
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II. Bewegungsgleiehung 
Die ‚Strömung wird durch die Bewegungsgleichung 


2 ı\2 
a BSG, ZEN) a But St 


i (Br ya ayl \ayy) 
und die Kontinuitätsgleichung 
» ou MW 
ar a 


beschrieben. 
Transformieren wir auf die dimensionslose Variable 


x 2 ect BR 0 SLR 
j Pe OR ö’ (), ne &). FL: 7). (= 3) 
(‘ Ableitung nach x) 


und eliminieren » mittels der Kontinuitätsgleichung 


# 


Ya E 
> 

| 

Bi 


e so erhalten wir — mit = er 
| L TE u 
120,0 l 2. ale 27 22 | 
in) 
va Fanlla ) = ai 3 Wi Zn B% 
w wo | Ro, — z 
v 


Infegreiiion über , wobei partielle Integration zur Anwendung kommt, ergibt 


Die 'Integrationskonstane 6, (x) hat dabei ende Bedeutung 


ö Ben, L Schubspannung an der Wand). 


; $ Vorausgesetzt, daß wir > die exakte, für n—0 ehnhdends; Lösung der Differential- 
EN gleichung vorliegen hätten, so würden in derselben zwei Unbekannte auftreten : die Integrations- 

i onstante ( (bzw. c;) und die Grenzschichtdicke 6. Zu ihrer Bestimmung würden dann die beiden 
4 n en, die wir aus den Randbedingungen € 


. (9) 


gi: Bi 


’ v zen fen) Pan and Sea 


+) Ge) ea Hi fen-| wu) ” [oin N : Ä O 


ADS TERR a le) 2a EN 
Re Re\on), ow\y), 


T L 
x _ Definiert man als, örtliche Widerstandzahl = ; °., so haben wir 
Fi ER Ä 2 RES 5 a , j y 
Re Ber N 
5 sr * Re wa n ur, _ aÜ= ) a 


ER; 
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Diese Gleichung spricht den Impulssatz für die ebene Platte der Breite Eins aus. Dieser besagt 
für die ebene Platte, daß der aan" gleich dem Impulsverlust ist: 


W(e) = Teutm—n)ay 
(z von der vorderen Plattenkante aus en) 


% r x 
Mit W(a)= [t,dx (wir betrachten nur die eine Seite der Platte) ergibt das 
0 / 


d ö(z) 
| ut uay 
bzw. 
d ö(x) 
PAIR ehe 
PR TE a “av. also (10). 
0 


Die Randbedingung o=1für n„=1 und der Impulssatz sind offenbar äquivalent dem Gleichungs- 
system (9). In den Abschn. IV und V werden wir, indem wir die genannte Randbedingung und 
.den Impulssatz auf die uns zur Verfügung stehende Näherungsfunktion für die Geschwindigkeits- 
verteilung anwenden, die Grenzschichtdicke ö und die Widerstandszahlen c; und c, näherungs- ° 
weise bestimmen. & 


III. Integration über der Reibungszone der Grenzschicht 
Bei dem jetzt erfolgenden ersten Integrationsschritt (s. Einleitung) verweisen wir auf den 
analogen Schritt bei der Behandlung der Rohrströmung [1]. 


In der Differentialgleichung (6) können wir, bei hinreichend kleinem n, gegenüber der 
Integrationskonstanten c, (6, #0 vorausgesetzt) die ee 


1 o [on [on}) +32 [oan 


vernachlässigen, wenn man beachtet, daß an der Wand wegen der Haftbedingung —=0 und 
9%w 


Fr 


0) 


Öx 


— 0. (Eine Abschätzung der Größenordnung der Trägheitsglieder erfolgt in Abschn. IV, 4.) 
Weiterhin beschränken wir uns auf ein Gebiet, für das 

l | | | 

5 ln Mm) Ders ia 2 MAR DTRER e ONBEeAe (11) - } 

(7. Dimensionslose Dicke der laminaren Unterschicht). 


Für den so definierten Bereich der n reduziert sich dann die Differentialgleichung (6) auf 


1 dw neh f 
FR un ar EEE 12); 
Re; on a Ö on a ( ) 
Wir erhalten dann in analoger Weise wie in [1]: 
ah u f 
1% 2 BERN. Er. 
a Y1+4x Rec, (n no)? ) 


„a (13) 
2 1n {2x Ro Yo (nn) + VIER RE cn — —} + Recon Ber 


für 0< n—m<l, d.h. ir * (n— n,) und Trägheitsglieder gegenüber c, vernachlässigbar. 
une 


Wir transformieren auf die Schubspannungsgeschwindigkeit v, = 7 | x A 
u 4 
VE U Ya 


und 


“ ne ee A re PER te Pc Fr 
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ı 


Die Geschwindigkeitsverteilung lautet dann — mit = %Y und yv-, == 
4 %% 
11—-Y1+42: x — 
ee u an a TFRREE + (16) 


ie 2 (Xx—X0) 
für 


0<y—x<Reaya. 
Für die laminare Unterschicht selbst haben wir 


VE un Re N a EEE Sen (17). 
Für größere n erhalten wir wieder die bekannte Gestalt des logarithmischen Verteilungsgesetzes: 


om (Refen) + a [14 mAr+rRajanı] | 


für ae ERIDLR IE 
- 1 
—ı<£n<i 
Re | 
bzw. 
1 1 
vorn tin n®x 8X] 
für nn RER Br TER A BE. 
} 5, <xx<Reo ya 


Dieses logarithmische Verteilungsgesetz, das zuerst von L. Prandt] [2] abgeleitet wurde, findet 
sich für im Bereich der turbulenten Reibung gelegene wandnahe Punkte experimentell bestätigt 
(s. Bild 1). Schultz-Grunow [3] gibt an: 


; EB TATEN ng ea Aa at (20). 
Für unsere Konstanten ergibt das 
= gg. ET EEE TRNE (21) 
Bes 
2, t,mia=40, 0 md. 2 


Für große Reynolds- 
sche Zahlen Re, = Hı® 


entspricht der Grenz- 
schichtdicke ö6 nach *# 
Schultz-Grunow [3] 
ein x, für das Igyx bei 
3,öliegt (s.Bild1). Man 2 


) sieht, daß dem gegen- 
über x, eine sehr kleine 

Zahl ist. “ 

Zusammenfassend 

können wir ganz analog en 

der Rohrströmung fest- < Schultz- Grunow, Lufo 17 (1940) 

stellen: Unmittelbaran ” laminar 

der Wand haben wir Theorie (lg Rey = 6,8) 

F laminare Strömung mit EP an een logarithm. Verteilungsgesetz 

pie dem Verteilungsgesetz 5 für wandnahe Punkte 


(17). Die dimensions- 
lose Dicke x, der lami- 


Er naren Unterschicht be- 4 
915,515. 0 05. 1,0 15 2,0 2,5 3,0 (2) 4,0 
N } j Es folgt ein Be- Bild 1. Geschwindigkeitsverteilung an der ebenen Platte 


reich, in dem laminare 
_ und turbulente Reibung von gleicher Größenordnung sind, mit der Geschwindigkeitsverteilung 
_ nach Formel (13) bzw. (16). 
Der obige „‚Übergangsbereich“ führt in den Bereich der turbulenten Reibung der Grenz- 
TER schicht. Es ist bekannt, daß dieser Bereich bis dicht an die Wand reicht. Mathematisch haben 


’ ar 
EM 

er iur. 
er 


kun 
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wir.für seine untere Begrenzung 7 > 2 gefunden, was mit 5 {= 3,5 für die Grenz- 
ER 


schichtdicke bzw. nach (15) mit Res Yc, = 3160 und x — 0,388 ie Korn 
ergibt sich nach Schultz-Grunow [3] (s. Bild 1), daß die untere Grenze wie bei der Rohr- 


ergibt. Experimentell 


strömung bei a, = 50 gelegen ist, was n == 0,016 entspricht. Die Integration erstreckte sich 


“ in diesen Bereich hinein, wobei eine obere Begrenzung gegeben war durch die Gültigkeit des 
linearen Ansatzes 1 =xy für die Mischungswegverteilung und durch die Vernachlässigbarkeit der 
Trägheitsglieder gegenüber der turbulenten Reibung. Die Geschwindigkeitsverteilung in diesem 
Bereich wird beschrieben durch (18) bzw. (19). Nach (18) hat die Verteilungsfunktion die Gestalt 


ya 
W= mn nt const. 


Man kann dieses logarithmische Verteilungsgesetz, das wir aus der allgemeinen Formel (13) durch 


1 
die Einschränkung n7 > TEEN gefunden haben, auch sofort aus der Differentialgleichung fi 2) 
x he; YCı 


gewinnen, wenn wir in dieser die Jaminare Reibung gegenüber der turbulenten Reibung vernach- 


I 
lässigen und den Ansatz zen machen. Die Berücksichtigung der Trägheitswirkungen ergibt, _ 


wie wir im nächsten Abschnitt sehen werden, ein in nlineares Glied, das für kleine »] von kleinerer 
Größenordnung als das logarithmische Glied In n) ist. Inwieweit das logarithmische Verteilungs- 
gesetz (19) die Geschwindigkeitsverteilung über der gesamten Grenzschichtdicke approximiert, 
kann aus Bild I nach Schultz-Grunow [3] ersehen werden. 


IV. Fortsetzung der Integration 
Wir gehen nun daran, die Verteilungsfunktion (18), die uns — wie oben bemerkt — bereits 
in den Bereich der überwiegend turbulenten Reibung führte, in einem zweiten Integrations- 
schritt (s. Einleitung) weiter in diesen hinein fortzusetzen. 


Wir haben es hier, bei Vernachlässigung der Jaminaren Reibung, en: (6) mit der Gleichung 


Ba le [en [u)l-siefen] 


zu tun. Integration ergibt 


n > 
lem Ya Fr: [ ] 

a ı dn . Ge 2 ee (23), 
wo | ] die obigen Trägheitsglieder bezeichnet. Wir erweitern den Integrationsbereich nun 
insofern, .als wir 

gen HUB. TA RER EEE (24) 
setzen und die Trägheitsglieder berücksichtigen. Wir wollen auch noch folgende Lincarisierungen 
zulassen: 

Ya —[]r Ja (i=3 5 N [ ) F (c, + 0 vorausgesetzt) 
£ 1 i j - ‘2 
A Da ı | 
3 ’ (1 Er E n) . 
Bin „ntu .n r 


Es folgt dann aus (23) 
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Der Anschluß an die beim ersten Integrationsschritt erhaltene Funktion (18) erfolgt in der 
Weise, daß die obige Integrationskonstante der in (18) auftretenden Konstanten gleichgesetzt 
wird! Das gibt 5 


% EEE 
const — 1 In (Re yo) + U 1—14Indx tz] N IE NTORN 
Es treten also in der Fortsetzung die Glieder 
j n ; 
BE 2 MO 
x | ES: a)JN ga | 
- 0 ü 
hinzu, die für kleine 7 von kleinerer Größenordnung sind als Inn. Wir beschränken uns 


im folgenden auf die Berechnung des Integrals | am, das für kleine 7-Werte von höherer 


N 
Größenordnung als das Integral = | []dn ist. 


ö 
Wir ermitteln nunmehr die Geschwindigkeitsverteilung innerhalb des gekennzeichneten 
Bereiches durch einen Iterationsprozeß. Dabei kommen eine Randbedingung und der Impulssatz 
zur Anwendung!). Infolgedessen wird unser Ergebnis, da wir es nicht mit der exakten Lösung 
) zu tun haben, es einem gewissen Fehler behaftet Ei 


1. Erste Iteration: Wir gehen i in das tat | — din, das die Trägheitswirkungen wieder- 


gibt, mit der Funktion (18) ©, = —- Inn + const a Funktion (18) stellt ja die Geschwindig- 


keitsverteilung in einem A eatück des Bereiches der turbulenten Reibung dar (Abschn. III)?). 
Man erhält dann, wie man sofort übersieht, Funktionen von n mit Koeffizienten, die — außer 
von x — von c, und ö abhängig sind. Zur Berechnung der Funktionen von c, und ö ziehen wir 
(vgl. den Hinweis Abschn. II) die Randbedingung —1 für „—=1 sowie den Impulssatz heran. 

Unterwerfen wir unsere Funktion (18) der Randbedingung &,—=1 für n=1, so erhalten 
wir die Gleichung 


ı- Um (Re, ya) a TA oe 


Diese Relation ist natürlich mit einem gewissen Fehler behaftet, da (18) ja nur in einem gewissen. 


. wandnahen Bereich die exakte Lösung darstellt. Wir können für (18) dann auch setzen 
[z 
ot nn A OT ED (28). 


Wir führen die Berechnung mit dieser Funktion durch und erhalten 


2 ; [tax ae) +('*) 38) - 301% en 
| era bi rl) men (2 


a = 
#4 (' Ableitung nach «). J 


H); Man kann den folgenden Iterationsprozeß mathematisch so beschreiben: Ausgehend von der 


\ ' Funktion (18) RL In 7 + const wird die Funktionalgleichung (25) 
BT}, : n 

Ba, Ve, u Im} 

a Dort : En Ense rg FR er - const 


 iterativ unter Beachtung der Randbedingung ® = 1 für n = 1 sowie des Impulssatzes gelöst. Dabei werden 
Glieder die für n<1 von kleinerer Größenordnung sind als n, vernachlässigt. Diese Vernachlässigung ist 


’ 
geboten. durch die Beschränkung auf den Ansatz (24) 5 = x + „n? für den turbulenten Mischungsvorgang. 


Mn =: Für die untere Grenze des Bereiches der BR Reibung hatten wir in Abschn. III die Bedingung 


n> ie Diese Zahl konvergiert für Res — oo gegen Null, wie aus (27) folgt. 
a Fer a 


DE San a are. Se) 20 


PESSE 


TI 


- 
a 


a 


RR rt ER 


| an "erhalten wir C, als universelle a * nicht von x abhängt: 


a A 
£ 
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Nach (25) gibt das eine Entwicklung 


; a (++ Buononlan—Qm+teonst 22... (29), 
wo ‘ 
a EN! al) 2 7 RER (30), 
B= An TR II RR er See (31). 
ar (alet ey erstere) am 


Die Berechnung der Koeffizienten A,, B,, C\ erfolgt so: 
%) @—=1 für n=1 ergibt nach (27): 


lea) + EI Hnduhez Eu, 22 A RE RN 


Diese Formel liefert eine Aussage über die Abhängigkeit des örtlichen Widerstandskoeffizienten 


% (4-54) von der Reynoldsschen Zahl Re. Es folgt aus ihr, daß 


c, = 0 für Ra a \ 
Für große Reynoldssche Zahlen Können wir also c, als gegenüber 1 kleine Größe behandeln. f he 
Durch Differenzieren — x, wird als Konstante behandelt (vgl. Einleitung) — erhalten wir 5 
aus (33) 
0 ”. N 3 
oe) N als a) SENT} se RE (35) ” 
% % % 33 
und oe 
Ya Ya Yan 
er IE 1181. | 12.2, CE Sa Are ae 36). 
eine U ee IK (86) | 
ß) Der Impulssatz (10) ergibt mit (28) und unter Berücksichtigung von (35) : A 
—\, = | . 
Ya Ya, 098 We 
«=(0Te) a2 ee Re viazk 


y) Für (, folgt dann aus (32) mit Hilfe von (86) und (37) u nr Eu 5 


C = n a 3 . 
Für A, nach (30) und mit Hilfe von (35) und (37) 


Entsprechend 


gegen Null Veaenibaner wir die für große Re a 
sehr kleinen Glieder mit den Koeffizienten 2 und B,, so ‚liefert die c BI 
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Das zusätzliche Glied —C;n gibt den Einfluß der turbulenten Rejbung in zweiter Ordnung 
gemäß dem erweiterten Ansatz (25) = =xnH+ un? wieder, sowie den Einfluß der Trägheits- 


glieder i in erster Ordnung. Beide miteinander in Konkurrenz stehenden Einflüsse liefern das 
. Inn lineare Zusatzglied, das für kleine 7 gegenüber dem logarithmischen Glied von kleinerer 
Größenordnung ist. 
2. Zweite und folgende Iterationen: Die zweite Iteration verläuft ganz analog. Wir 
beschränken uns daher auf die Wiedergabe der Ergebnisse: Unterwerfen wir die Funktion (41) 
der Randbedingung w, = 1 3 n= 1, so können wir auch setzen 


Br rele  {mn+G(d—m)}. 
. Mit dieser Funktion erhalten wir dann nach (25) eine Entwicklung 
2 3 
la hanlır 4n4 Bund DE .|- — 0m +05 +H% |+oonst . . (42). 


Vernachlässigen wir die Glieder, die für wachsende Reynoldssche Zahlen gegen Null kon- 
vergieren, so liefert die zweite Iteration 


1a TInn—Qn+@ 14 const 


y A PER 2| 
mit 
A LEN, A 
ns nen er 


2 


Das in n quadratische Glied ist für kleine n von kleinerer Größenordnung als das in n lineare 
Glied und wird vernachlässigt. Es möge bemerkt werden, daß, wie man leicht aus den Entwick- 
lungen auf S. 314 sieht, die Erörterung eines in 7 quadratischen Gliedes eine Erweiterung des 
Mischungswegansatzes auf drei Glieder erfordern würde. Wir erhalten dann 


c 1 
| = Inn — O;n} + const, G=.+ DT ae (43). 
Be Man sieht nun schon, was weitere Iterationen ergeben. Wir erhalten 
2 Ahr. {In n— Din VE const so na ei (dA 
| mit 
. WETTER, NHL 1 2 ER 44 
te Beta (44) 
a Es folgt 
277 UBRAHE 
ns £ > TE R En On * 120° 
was 
| 1a re 45 
0=1+5 + = N ee 2 1 RR (45) 


er; gibt. ‚Im Endergebnis erhalten wir also in Fortsetzung der beim ersten Integrationsschritt 
we Funktion (18) 


ı 
re. 
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Insofern (' eine universelle Konstante darstellt, besitzt die Funktion universellen Charakter, 
d.h. sie ist kennzahlunabhängig und variiert nicht von Strömungsquerschnitt zu Strömungs- 
querschnitt. Ein solches Verhalten zeigt nun auch tatsächlich die experimentell gefundene 


Verteilungsfunktion “4% gn der ebenen Platte (s. Bild 2 nach Schultz-Grunow [3]). 
Ox 
3. Bestimmung von u/x: Bei Transformation auf die Schubspannungsgeschwindigkeit v, 
erhalten wir aus (46) unter Benutzung der Transformationsformeln (14) und (15) 


C 1 
y-llng- > E bonst Ts EN Te TERRA]: 


In der Fortsetzung verliert also Formel (19) durch das ortsabhängige Zusatzglied ee X 
ö 
ihren universellen Cha- 
rakter. Das zusätzlich 
erhaltene Glied kann 
geometrisch so inter- 
pretiert werden, daß es 
bei der logarithmischen 
Auftragung der Ge- 
schwindigkeitsvertei- 
lung die Abweichung 
der Tangente von der 
Geraden, die in dieser 
Darstellung die loga- 
rithmische Verteilung 
repräsentiert, wieder- 
geben muß. Das gilt 


Experiment auch für (48).. Anhand 
Schultz - Grunow, Lufo 17 (1940) : j derin Bild 1 und 2 dar- 
en gestellten experimen- 


tellen Geschwindig- 
keitsverteilung haben 
wir 


0 0,5 10 15 2+ 195 2,0 
Bild 2. Geschwindigkeitsverteilung an der ebenen Platte Ce 1r (50) 


bestimmt. In Bild 2 haben wir die mit C—=—1 berechnete theoretische Kurve (48) mit der 
von Schultz-Grunow [3] gemessenen Geschwindigkeitsverteilung verglichen: Der Tan- 
gentenverlauf wird danach für kleine 7 durch das gewählte C gut wiedergegeben. Die bei 
größeren n-Werten auftretenden Abweichungen halten sich in mäßigen Grenzen. Koordinaten 
der berechneten Funktion: 


n | 2+lgn | bc Eee n | 2+lgn ra a 
DJ Vox 
0,032 0,50 11,38 0,178 1,20 8,565 - 
0,056 0,75 9,85 0,316 1,50 4,725 
0,100 1,00 8,25 0,562 1,75 2,61 
In Bild 1 haben wir des weiteren die sich mit CO = — 1 ergebende Geschwindigkeitsver- 
teilung (49) für lg Re, — 6,8 mit lg e, = 7,5 — 10 aufgetragen. 
Re,Ya= Res) — 
. | Nach Formel (68)  #e&jcı” Re, Ba j 
Aus C=—.1 folgt nach (46): 
BI BB re (51). 


Eine weitere Vergleichsmöglichkeit bietet die Mischungsverteilung. Experimentell hat für diese 


g Schultz-Grunow [3] den in Bild 3 angegebenen Verlauf ermittelt. Mit »—0,388, — N 1,33 


© 


Ds N 
“ En . ’ * 

Er eh 4 £ -" ‘=, & 
a Gab y E s n ” Di he. hr 
5) n, ö ? 5 5 
. ® Nr m a N . R Hr” N 
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“erhalten wir aus (24) für: n=0,1 .— — 0,0336 (0,0388), 
n=02: = — 0,0570 (0,0776). F 
“ In Klammern haben wir die Werte angegeben, die der Ansatz 0.3 xn ergibt. 


ö 


4. Größenordnung der Trägheitsglieder in der Reibungszone: Wir können nun- 
mehr auch die Größenordnung der beim ersten Integrationsschritt (Abschn. III) vernachlässigten 
Trägheitsglieder berechnen. 

Mit (47) erhalten wir für die rechte Seite der h 
Gl. (6)°) für große Re-Zahlen /6 


ya ’ c Ye, Be 20 
(1) (n-5)+ (2) {nInn+Cn (i—n)} 
ann) nun a 


0,04 


) ra 11-—— nach (56b). era. Sa E 


Je Bild 3. Mischungswegverteilung an der ebenen Platte 2 
2 A nach Schultz-Grunow, Lufo 17 (1940) "u 


Die von n abhängigen Glieder stellen den Einfluß der Trägheitsglieder dar. Die Vernachlässigung 
der Trägheitsglieder beim ersten Integrationsschritt ist also gekennzeichnet durch die Bedingung 


Bl 
1 RO 12 
2 


N 
5. Verdrängungsdicke: Aus der Definition ö* — ö / (1 — w)dn folgt mit (47) unter 
ö 
Beachtung von (8) 


5 % 
Setzen wir ein nach (21): x—=0,388, nach (50): C = — 1, so Re — 2,73 Ye5, während Schultz- 
4 


Grunow [3] aus den Messungen 2, 36 Ya erhält. 4 


Bi), 6. Querkomponente: Wir berechnen v aus der Kontinuitätsgleichung (4) der Bewegung. 
F: Mit (47) erhält man unter Berücksichtigung von (35) und (36): 


BE lnnenor)] 
SER Seal (56b) | er RN | 
Der 


für große Re-Zahlen. Damit gewinnen wir 4 


v3 ar 1 B; 
PR 2 ig ao lin 3) u 


D) Bezüglich der Anwendung von (47) bis zur Wand n=0 hin ist zu bemerken: Aus der Anmerkung 
fS. 315 eh, daß die Dicke des Übergangsbereiches einschließlich der laminaren Unterschicht für Res ® 


— Laer Null konvergiert. Auch der Koeffizient Ye, in (47) konvergiert für Reg> x gegen Null, 


er; ; denn nach (27) oder (55) haben wir Ya 0 wie ER, Der dadurch begangene Feh- ? # A 
Er es Ycı j 
in den inerer Größe von Gl. (6) Ban zur Wand Ar =() zur Anwendung bringen, ist also für 3 3 
nn ö ne j' 
f i 22* nf 
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für große Re-Zahlen. Für den Betrag der Querkomponente am Rande der Grenzschicht (n=1) 


a: folgt ; N 
E | 2 % ; (53) 
m De 9 f . . . . . . . . . . . . 
E“ Mit C=—1 und bei Bezugnahme auf die Querkomponente v, am Rande der Grenzschicht 
Br erhält man i | 
Be - | 3 RR 
= »-3(n+®) n 
AR In Bild 4 und 5 haben wir experimentelle und theoretische Ergebnisse bezüglich der Quer- h 

% komponente miteinander verglichen. 

ah u a in ie in nn 25:10* 


Querkomponente am Rande 
der Grenzschicht 


v 
Verteilung " ur 


Re > 


15:70°* 


10:10°* 


. Messung Schultz - Grunow 
Lufo 17 (1940) 


Messung 


4 Schultz - Grunow, Lufo 17 3:9” h 
ya (1940) Theorie 
Theorie 


Pu 


3 
0 0.2 0,4 a: 0 6,0 6,2 6,4 6,6 IgRe, 70 
Bild 4. Querkomponente der turbulenten Strömung Bild 5. Querkomponente der turbulenten Strömung 

an der ebenen Platte an der ebenen Platte i 


7. Wir geben vergleichsweise die empirischen Koeffizienten, die wir bei der Behandlung 
der Rohrströmung [1] erhalten haben, an: 


Rohr Ebene Platte 


Xo —r 7,01 Xo = 5,515 3 
”“ = 0,417 20,388 
ER di ee es 
% x & 
AR V. Grenzschiehtdieke und Widerstandsformeln 


Wir berechnen die Grenzschichtdicke ö und die Widerstandszahlen c; sowie c;, indem wir 
die Randbedingung ®=1 für n—=1 und den Impulssatz auf die gewonnene Näherungsfunktion 
(46) bzw. (47) für die Geschwindigkeitsverteilung anwenden (vgl. den Hinweis in Abschn. II). 

Unterwerfen wir (46) der Randbedingung =1 für n=1, so R 


. 


} 1 1 1 ; - u“, ah 1% 
nn a a oe ER N. 
M ' ’ BL 
Der Impulssatz (10) ergibt mit (47): / H TR SNSDE en 


wobei 


} + Zu a 
3 Es Be - 


AR ex 
Bo Be tath. Meöh, 1961 Kork r 1 ewski, Becher der erhalten BES längs der Fee Platte en 


BR 1. Grenzschichtdicke : 
Aus (55) folgt sofort a 
R,=- an, 9 ld AR MER) K 
en. ya, Kr 
| ” f 
womit man die dimensionslose Dicke Re, als Funktion von c, gewonnen hätte. Bin 
Setzt man (57) in (56a) ein, so 0 
£ d _ % Ya 
& «Ve NL N ee ent 0 RER OR a 
= 4 dRe, 1 ii BE Four? ) (58) Be 


mit der Abkürzung 


121 e 
N TE 1+0—x De. 
le 4 2 eh: hr “FR 


Integration ergibt, wenn wir die Forderung ö—>0 für 2—0 stellen, > 
; » | A „Va h 1a Ri: 
ER = — = a | 61 B 
2 | Re, Br a EN ON een m 
) i MER 3 i 
vi A 
Damit haben wir auch den dimensionslosen Abstand Re, als Funktion von c, gewonnen. \ 
Gl. (57) und (59) stellen eine Parameterdarstellung der Funktion ö(x) dar. Wir merken Br 
noch für späteren Gebrauch an F 
- Er 
2 .e 
FB 1. ea ie ee (OP 
C & er a 
ger i 
f In der Asymptote (B< < 1) läßt sich das Wachstum der ERERCHTANETICKE ö mit zunehmendem # 
. ' Abstand x von der Plattenvorderkante explizit beschreiben. ” ; 
3 £ Bu: i Bi 
Be j Ae’a . Br. 
2 N; | Nach (59) by ‚Re, ar etc hie Et (59a). R: 
4 N 6 r Re; 2 Ya 60 | Er 
ach (60) Dabei BEN; RN N a). 
Ana en folgt: In Re, = Z j vi damit aus (60a): i 
vr 
„2 Re, De; 
Re m 6 In Re, NL N EIER (61 a) r j: 
„2 & 2 
Bu... . (61b). er 
e Er c In Re, $ RR Be: 
2 ö Sad 
zn un. 4 ‘ ® R = y 
Ilse in der Asymptote wie während man bisher bei der ebenen Platte ein Wachstum : 
en zschichtdicke ö mit 2*/5 angenommen hatte. | E 
j 2. Widerstandszahlen A 
ö 
ee ER 1 (62). X 


; FA 
Pro. 
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wo W einseitiger Widerstand einer ebenen Platte der Breite Eins Lad der Länge . 


= Aus 
$ b % % s 
BEN W=/[1da=ou:[cd« 2 
folgt mittels (56a) 
Bar.: 22a + le SR a En (63) 
Er Setzen wir für ö/& Formel (60) ein, so erhalten wir bis auf Glieder kleinerer Größenordnung 
% F ‚ rn 2 Cı 
Br - ‚= ! 
e 1“ +,% & j 
ae, a»? 
ke | Die Auflösung dieser Gleichung nach ec, bzw. 4—=2c, ergibt unter Benutzung des Zwischen- 
x ergebnisses 
r 1. 
N EI 9 MR, ER 
Bo: m y2 rzE as = ae (64) 
Ber, She 3 
an (hierbei wird 2 6; # x? benutzt, ein mittlerer Wert von c; ist 0,002) 0) 
28 1 ; 
Bi . j E Re \ y Bu . 
, d G — 1 G; y,-+ Pe N Beh (65). { vr 
- Damit haben wir eine Aussage über den N zwischen der örtlichen Widerstands-- r 
zahl c; und der Gesamtwiderstandszahl 10 gewonnen. Im Vergleich dazu gibt Schönherr [4] 4 


die Formel an 


während Schultz-Grunow [3] 


mit 
0,427 
= 0,407 + 1g Re, ' 
Der Vergleich der theoretischen Formel mit den experimentellen Formeln ergibt 
für Re, = 10° mit c;— 0,00294 für Ro = 10° mit c, = 0,00145 
nach Schönherr *7 c; = 0,00246 c = 0,00127, 
nach Schultz-Grunow c = 0,00242 = 0,00126, 
Formel (65) c; = 0,00239 Aha! 4 = 0,00125, 


3. Widerstandsgesetz all: Br 
Wir ermitteln nunmehr das Widerstandsgesetz, d.h. die Abhängigkeit der Gesamtwider- | io 


standszahl von der Reynoldsschen Zahl a Ars Ban (65) haben wir he 


in (Ralo)tz 2 14mar+ 0). 


Wir werden die gewünschte Formel erhalten, wenn es uns gelingt, a und Res durch or und} Be,. 
zu ersetzen. re nr 


Br 


Bi a) Zunächst haben wir 2 
N: je Roya = Dr A 


’ was mit (63) ergibt 


+. 


r & d But 
KR er A "7. , n 
a #. a N ARM EN = Bi a Fu RUE 


Dir bei Aa 


A Bein nn Okt. Br Szablews ki 1% Berechnung Br leuten ‘Strömung längs der ebenen. Platte 323 
I 


| ß) Nach (65): 


Ya: Ye y2 W; 
\: ot 3 530 
Aus der Entwicklung in eine unendliche Reihe erhalten wir 
a) el | 1 
lee bwin — Bun I Re, 69). 
y) Ferner haben wir Be, 
In‘ + Fa, FA lnard, ar ae! BE: 1 . 
a; 2 M x &) 2 f“ \E 
Ersetzen wir hier nach (64) ia durch + Va- E T, so hi 
IK a“ RATE 2 | RÜ\ 
) Zr a | Be 
G vaw-|er Ag 0% %\11 2 Er 
Be Int PP un AT ya aa lat 2. (70). E 
Unter Benutzung von (68), (69) und (70) ergibt sich dann aus (55) g 1 2 
ERS, LE, & | | Br: 
1, 21 ya le in 2 
= SE Ru)+,. 2+ ma) -C+rn al 0 


„2 5 | 


. wobei nach (56): ; 
NER = — 2430-50. 


2 ’ 
A 1 ; - 
| Nenlabhlässigen wir noch auf der linken Seite der Gleichung die gegenüber — sehr kleinen . 
übrigen Glieder (vgl. die unten angegebenen numerischen Werte der Koeffizienten), so erhalten A. 


wir schließlich 


82; | HL Ber » 2 Fe en 1 


j 22 £ 2,3 > 
Ele" re. Sn Er are Ten TEEN garten, :*; ( Si 
Br a lg (Re, +2 +m@M) O4 rn — in a 2. 

en | 2,3 a 2 ee 

Die Koeffizienten werden bestimmt durch die dimensionslose Dicke x, der laminaren Unter- SE 
schicht und die Koeffizienten x, w der Mischungswegverteilung. Wir hatten gefunden in (21) 72 
und: ee 7 
F CR ie : = 0,388, Din RI 


Be De N 
x 


= — Il en ee nen. (2) K 
PER Y9=g(Re, ) — 0,202 en) | 
gleich führen wir die Widerstandsformel von Schönherr [4] » 
RE) 242 
a ler lg (Re, 6). R 
Schultz- Grunow B] angegchene Formel ’ K 
Baader nz | n 
ER 0, 407 a I SR a us Be. 
BR ; 
23 
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Der Vergleich der theoretischen Formel mit den experimentellen Formeln ergibt 
für Re, = 107 für Re, = 10° 


nach Schönherr c;— 0,00293 c; = 0,00156, 
nach Schultz-Grunow c, = 0,00294 ce, = 0,00145, 
Formel (72) c; = 0,00308 c; = 0,00160. , 


VI. Zusammmenfassung » 


Unter Zugrundelegung des Prandtlschen Schubspannungsansatzes werden die Differential- 
gleichungen der turbulenten Grenzschicht an der längsangeströmten ebenen Platte in zwei 
Schritten von der Wand her integriert. Für den übrigen Bereich der Grenzschicht wird die 
erhaltene Lösung als Approximation herangezogen. Es wird die Geschwindigkeitsverteilung für 
die Längs- und Querkomponente berechnet und eine Widerstandsformel aufgestellt. Ferner 
wird das Wachstum der Grenzschicht in der Asymptote berechnet. Bei passender Wahl der 
empirischen Koeffizienten der Theorie erhält man gute Übereinstimmung mit den experimentellen 
Ergebnissen. 
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KLEINE MITTEILUNGEN 


Erzwungene Schwingungen des linearen 
Schwingers bei nichtharmonischer Erregung. 


Bei der Behandlung der erzwungenen Schwingungen 
wird in den meisten Lehrbüchern nur der Fall der 
harmonischen Erregerkraft ausführlich behandelt und 
der allgemeinere Fall einer periodischen, nichtharmo- 
nischen Erregung mit dem Hinweis auf das Fourier- 
sche Theorem und das beilinearen Systemen geltende 
Superpositionsprinzip abgetan. Damit wird zwar das 
Problem im Prinzip gelöst, jedoch ist die zahlenmäßige 
Ausrechnung meist mühselig und zeitraubend. In 
Sonderfällen kann man jedoch auf völlig elementarem 
Wege zu Lösungen kommen, die weniger Aufwand zu 
ihrer zahlenmäßigen Auswertung erfordern. 

Wir betrachten den linearen gedämpften Schwinger, 
dessen freie Schwingungen durch die Gleichung 


ai ms +-ds+ c:=0.... (]) 


beschrieben werden. Durch Einführung der neuen 
Zeit 


pl 


6 
r=l\-I 
m 
läßt sich (1) in die dimensionslose Form 
dx ds 
APR + 2D Fr + EN) Eee (2) 
bringen, wobei B 
ee 
2 Yme 


eine dimensionslose Größe (das Lehrsche Dämpfungs- 
maß) ist. D=0 gibt den ungedämpften Schwinger, 
D=1 den aperiodischen Grenzfall. Die Lösung von (2) 
lautet: 

s=nac ?roos YL—Diit—r) .....08). 


Die Bedeutung der Integrationskonstanten v4 und r, 
geht aus Bild 1 hervor. Die Zeiten, zu denen die 
Schwingungsamplitude Maxima erreicht, folgen aus: 


D 
. (4). 


tg 1 — D? (tman Da ne ER 


Die Maxima selbst sind: i 


nn 
—D (mar+ en 


en 
mas = al — D?e 1 


n=0, 1,2-+-(5). 


Diese bekannten Beziehungen über die freien Schwin- 
gungen sollen dazu verwendet werden, die Lösung 
auch für erzwungene Schwingungen aufzubauen. Es 
wird dabei — wie das üblicherweise auch bei harmo- 


nischer Erregung geschieht — nur der eingeschwungene 
Zustand betrachtet, der sich nach Abklingen von 
Eigenschwingungen im Laufe der Zeit stets einstellt. 


1. Erregung durch eine periodisch wechselnde 
Kraft konstanten Betrages. 

Wir wählen den Betrag der Kraft so, daß die durch 
sie bewirkte Verschiebung der Gleichgewichtslage des 
Schwingers gleich der Einheit wird. Die Zeit einer 
vollen Perle des Kraftwechsels seit. Dann en 
die Gleichgewichtslage 2,1 jeweils nach der Zeit 7y/2 


O 


3 
3 
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von —1l nach +1 oder umgekehrt (Bild 2). Die Auf- 
gabe besteht nun darin, in die mäanderförmige Kurve 
der Gleichgewichtslage solche Teilstücke der freien 
Schwingung von Bildl hereinzulegen, daß eine 
stetige, kniekfreie Kurve entsteht. In Bild 2 ist das 


Al. lr) 


% 
0<V1-0°r,<22 


ZRr<VI-DU Tr, <a 


Lore 


4x <V1-02 ,<6r 


Bild 2 


für drei verschiedene Werte von r, angedeutet. Be- 
zeichnet # die um den Betrag 1 verschobene Schwin- 
gung, so muß gelten: 


eg: @ 


2 3 4 
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Durch Einsetzen von (3) folgt aus (6) eine Bestim- 
mungsgleichung für 75: 


tg ylL — Du, = 


T, 
Dr? 
ab a en 
D+e (VTDisinyr= 32 | Deos vi=D:7) 


nm 


em: 2 T T 
NDR (ve cos /1—D? o. —Dsin Y1—D? 2) 


(8) 


und aus (7) eine Gleichung für 24: 
2 
au= EE9)3 
A m J (9) 
cos YL— Dry +e ? cos yl- (9-7) 


Damit sind die Kennwerte der Schwingung bestimmt. 
Zur Berechnung der Resonanzkurven kann man nun 
wie folgt vorgehen: Für gegebenes D und r, wird 
aus (8) 7, bestimmt, dann aus (9) 74, aus (4) Tyan und 
aus (5) Zmax: Schließlich wird max aus tman durch 
Hinzuzählen bzw. Abziehen der Einheit gebildet, je 
nachdem, in welchem Bereich 7, liegt bzw. welches 
der beiden ersten Maxima Tas den größten Betrag 
für max ergibt. Das Ergebnis der Auswertung zeigt 


Bild 3. Dort ist als Abszisse ö -,2 ‚ also das Ver- 


hältnis von Eigenfrequenz des Schwingers zur Erreger- 
frequenz aufgetragen. Das ist im vorliegenden Fall 
zweckmäßiger, als die übliche Auftragung der Reso- 
nanzkurven über dem Kehrwert 1/6, also dem Ver- 
hältnis von Erregerfrequenz zur Eigenfrequenz, weil 
so die Maxima besser getrennt werden. Zum Vergleich 
ist in Bild 4 das bekannte Resonanzgebirge für har- 
monische Erregung in der gleichen Weise dargestellt 
worden. Es zeigt sich, daß die Werte des ersten 
Maximums von Bild 3 rd. 20% größer sind, als die 
von Bild4. Für ö<1 ist zwischen beiden Resonanz- 
gebirgen kein wesentlicher Unterschied festzustellen. 
Dagegen sind die Unterschiede für ö >1 beachtlich. 
Sie bestehen erstens im Auftreten weiterer Maxima, 
die allerdings nur bei geringer Dämpfung merklich 
werden, sowie in der allgemeinen Hebung des Ampli- 
tudenniveaus, das auch im Grenzfall 6 — © über 1 


de —————— 


a Bild 3 


Pr er 


326 


bleibt. Das ist physikalisch sofort einzusehen, denn 
auch bei außerordentlich kleiner Frequenz des Kraft- 
wechsels wird das plötzliche Umschalten der Kraft- 
richtung stets zu einer Schwingung anregen, die über 
die Gleichgewichtslage hinausschießt. Bei sehr lang- 
samen Kraftänderungen — wie sie bei harmonischer 
Erregung kleiner Frequenz auftreten — kann dagegen 
die Amplitude den Wert 1 nicht übersteigen 

Die Lage und der Betrag der Resonanzmaxima läßt 
sich leicht errechnen. Ein Maximum tritt stets auf, 
wenn mit der Kraftriehtung auch die Schwingungs- 
richtung wechselt, denn dann wird — verglichen mit 
den möglichen Nachbarbewegungen — am meisten 
Energie in das System hereingepumpt. Aus (8) und (4) 
folgt, daß dies der Fall ist für 


an yı — D2 


Tp n > 


r 0 — 
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ferner aus (4): ma =utNT, 
und schließlich: 


n 2, 3 


1 


Tp 
cos — 
4 


74, 


mas =tı FT 1l= 


wobei das Minuszeichen nur für ö</1 gilt. Die ver- 
schiedenen Äste der Resonanzkurve sind also für 
D=0: völlig gleich; selbst der Anstieg des ersten 
Maximums entspricht dem aller anderen, ist jedoch 
um den Betrag 2 nach unten verschoben. Die Dicke 
der Resonanzmaxima nimmt in der D-Richtung 
rasch ab. 


Im Sonderfall D=1 (aperiodischer Grenzfall) hat 
man als Lösung von (2) zu setzen: 


»=21(L Herner. 


3 


In den Fällen n=2, 4, 6, ... ergeben sich Minima 
der Resonanzkurven, weil’dann ein Größtbetrag der 
ber einer Hinschwingung hereingepumpten Energie 
bei der Rückschwingung wieder abgegeben wird. Die 
Extrema selbst haben den Betrag: 


Dxn 
2e re, 
Vextr = Damen n=inZ ee 
1 
se Falle harmonischer Erregung liegt das Maximum 
ei 


al 
yı —2 D: 


und hat den Betrag: 


1 
% sg nn 7 
a nee 
Zu einfachen Ausdrücken führen weiterhin die Grenz- 
23 fälle D=0 und D=]1. Im ungedämpften Fall D=0 
x tr. folgt aus (8): 
ER aus (9): 74 = z 5 
4 Tp 
cos 


Fa 
en 


N nee 


7 


Bild 4 


Die Bedingungen (6) und (7) führen jetzt zu den 
Werten: 


ER, 
Eee u > 
a 
2 ;z 
il 1er 
ar ( 2) 
% 2. 
— BE ’ 
» 
l+e (1402) i 
und schließlich zu: i Er 
c—1 
TER OT 
Imaz = BAce —1. 


2. Erregung durch periodische Stöße 
wechselnder Richtung 


+ 
A d 


 Z. angew. Math. Mech. 
Bd.31 Nr.10 Okt. 1951 


Für ED, 


. w 3 D 
a & 4 > > h { 
k 4 Pe A ’ ß 
- - we F u 7 ß Fi 
A >. £ FR vw 
FERN a a 2 Be; 
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Geschwindigkeit um den Betrag 1 verändern. Wieder 


T 
erhält man durch Aneinanderheften entsprechender = ar 7 
Teilstürke der freien Schwingung (3) die resultierende l-+e cos 1 — D? 2 
; y ERNST EN NIE 4 a 
‚ 0</1-02 ,<2x ee He 


X(T) 


T 


AK-U- DISK 


47<V1-0° 7, <6r 


Bild 5 


erzwungene Grundschwingung. In Bild 5 sind drei 
derartige Fälle gezeichnet worden. Die Übergangs- 


bedingungen lauten jetzt: 


at 


Durch Einsetzen von (3) folgen daraus die Integra- 
tionskonstanten 24 und 7;: 


» 


w. Br 
e ° sin -Dy 


Mu= 
1 


E9:2 
1 Eee N 
Yı-D2| sin YI—-D! x,—e "nid Fn)|. 


Das Maximum der Amplitude wird — wie im vorher 
behandelten Fall — aus (5) ermittelt. Für hohe Fre- 
quenzen, bei denen der Größtwert der Amplitude in 
den Knickstellen selbst liegt, also für 


0< mas <Z 


gilt einfach: 


Imaz = (0) = zu c08 yı-D Ton 


Das Ergebnis der Auswertung zeigt Bild 6. Dabei ist 
der Ordinatenmaßstab gegenüber den Bildern 3 und 4 
verdoppelt worden. Der bemerkenswerteste Unter- 
schied gegenüber Bild 3 besteht darin, daß die Reso- 
nanzkurven jetzt für ö&=0 bereits mit positiver Stei- 
gung beginnen. Das Zustandekommen dieser Tat- 
sache läßt sich sofort aus einem Vergleich der entspre- 
chenden Schwingungsbilder Bild2 und 5, obere Kurven, 
erklären. j 


Der Grenzfall D=0 ergibt jetzt: 


ae 
ag? für d<1 
45 g 4 ’ 
el 
A funm0>>1 
Tp 
208° 


a 
b 
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Im Grenzfall D=1 hat man: Das Maximum liegt bei: 

wi Me] Tmas = re 
B | | 
ro Fi Nr und hat den Betrag; . { : 
4 Er 0 mas | 7 
e”, 2 Tmas = UA ce . q 
at h 4 
a, ö 
BE: 2) Es wird jedoch nur im Falle 0 < man < 7 wirklich Ar 
in Tp 2 T 2 

ie BI erreicht. Für tyar > gilt: = > . E 
H EUR 2 l } 
ir DAT 7, 2, j E 
% r i + E Zmaz = (0) = 74. 4 
Ber | 
» 


l Sl 


TE ee | 


1 


30 ES . | & h 

g ! 7 
2 AN) er, 2 E ;* 
0 i 


10-8 
re." 
Ay 
i - 
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3. Erregung durch periodische Stößein 
gleicher Richtung 


Wirken die Stöße stets nur in einer Richtung — 
was in der Praxis häufiger vorkommt, als der unter 2. 
behandelte Fall —, so gelten die Übergangsbedin- 
gungen: 

2(0) = 2(%p), 


dx de 
= — I — —l 
dt )r=0 dr (1=7,) 


Daraus folgt in üblicher Weise: 


Dr 


12e' Posi Di 
na em, 
6,2 sin /1 — D?r, 
Ra 
1 
I D 


VYı-D: [sin Yı-D® + e Fein Yı-D:(t ey 


Die Schwingung ist jetzt unsymmetrisch zur Gleich- 
gewichtslage 2—=0. Sie besitzt eine mittlere Abwei- 
chung x, von dieser und eine halbe Amplitude x4 
gegenüber der um x» verschobenen Mittellage. Es 
wird gesetzt: 


. 


1 
m = 5 (Tmazı + Zmax2), | 
E77) =; (tmazı — Zmazx2), | 


wobei die tmas wieder aus (5) folgen. Jedoch muß 


in (10) an Stelle amaz2 stets (0) eingesetzt werden, 


falls 
|2mas2 — zm| < |x(0) — ml. 


Die Auswertung zeigt Bild 7 für x4, Bild 8 für a 
Dabei ist gegenüber den bisherigen Bildern der Ab- 
szissenmaßstab verdoppelt worden. Während das 
Resonanzgebirge früher bei jedem ungeraden ganz- 
zahligen Wert von ö ein Maximum zeigte, tritt jetzt 
bei jedem ganzzahligen ö ein Maximum auf. Diese 
Tatsache läßt sich physikalisch durch den Fortfall der 
Bremswirkung des Gegenstoßes im Falle 2 leicht 
erklären. 


In Bild 8 ist die D-Achse nach hinten positiv ge- 
zeichnet worden, um das Gebirge übersichtlicher zu 
machen. 


Wieder ergeben die Grenzfälle einfache Lösungen. 
D=0 gibt: 


Tp 
men 
x 2 
NE I 
#7 
2 sin 2 
er 
Beate ir für Tmar2 < 1p; 
9% sin 2. | 
2 
1 Tp 
iZm—z7tgz 


für Tmaz2 > %p. 


REEL EP ER EEE, 
IRRE EN tr, 
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Ostaschkow, UdSSR. 


Knickung elliptischer Ringe. 


1. Fragestellung 

Im Zusammenhang mit Fragen der Tragfähigkeit 
von Streckenausbauten im Bergbau interessiert die 
Kenntnis der kleinsten Knicklast eines elliptischen 
Ringes unter einer speziellen Außenbelastung. Man 
benutzt in der Praxis nämlich neuerdings elliptisch 
geformte, geschlosssene Streckenbögen und versucht, 
das Achsenverhältnis der Ringe so zu dimensionieren, 
daß der Außendruck des Gebirges im unverformten 
Ring keine Biegespannungen hervorruft. Nimmt man 
an, daß dieser Idealfall tatsächlich verwirklicht ist, so 
ist damit, wie sogleich gezeigt wird, eine spezielle Be- 


lastung Z = — ——— vorgegeben. Hierin bedeutet Z 


Ro(®) 
die in Richtung der Ringnormalen wirkende Last je 
Längeneinheit des Umfangs, — $ die konstante Nor- 
malkraft ($ < 0, Z> 0, für Druck) im Ring und R,(9) 
den Krümmungsradius der unverformten Ellipse, ab- 
hängig von dem Winkel # der Normalen gegen eine 
feste Richtung (vgl. Bild 1). Die Berechnung der kri- 


Bild 1 


tischen Außenlast Z;läßtsich nach A.Lockschin[l] 
auf die Berechnung der kritischen Normalkraft $7 als 
Eigenwert eines Eigenwertproblemes mit gewissen 
Besonderheiten zurückführen, 


2. Herleitung der Beulgleichung 
In den Gleichgewichtsbedingungen für ein ebenes 
Element des gekrümmten Stabes mit fehlender Tan- 
gentialbelastung (vgl. z. B. [2]): 


9+Q9=0, 9! —-S-ZR=0, M’—-QOR=0..(l) 


werden außer der Querbelastung Z und dem Krüm- 
mungsradius R des verformten Stabes auch die Nor- 
malkraft S, die Querkraft Q und das Biegemoment M 
als Funktionen von ® betrachtet und die Differen- 
tiationen nach # werden durch einen Akzent be- 
zeichnet. Zunächst verlangen die Gleichungen (1) für 
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den biegungslosen Zustand aus der Forderung M(9) =0 
die schon genannte Normalbelastung Z= — RR. und 
0 
im übrigen ist auch Q—=0. Im verformten Zustand 
wird dagegen ein Biegemoment m(®), eine Quer- 
kraft g(®) und eine geänderte Krümmung A(®) 
auftreten. Mit der Grundgleichung der Biegelehre 
"(2 — = (EJ = const Biegesteifigkeit) folgt 
dann aus (1) nach Linearisierung für das Biege- 
moment m die Differentialgleichung 3. Ordnung mit 
periodischen Koeffizienten 


% Ave 3 % Ro 0 2 
(2) | =. Ir m) Er (2); 


In die Koeffizienten von (2) ist der Faktor « (kleine 
Halbachse der Ellipse) eingeführt worden, um auf 
dimensionslose Koeffizienten zu a Der eben- 
EG 
EJ 
wert, denn zu (2) kommen noch dreilineare, homogene 
Randbedingungen. 


falls dimensionslose Parameter } — ist ein Eigen- 


3. Randbedingungen 


Zunächst muß m(®) für einen geschlossenen Ring 
die Periode 2r besitzen. Es muß also für jedes 
gelten: 


mA—m(dA 27) „2. (2)’. 


Die noch fehlenden Randbedingungen werden durch 
die Periodizität der elastischen Verschiebungen der 
Stabachse geliefert. Wie in der Biegelehre für schwach 
gekrümmte Stäbe abgeleitet wird [2], kann man die 
Krümmungsänderung der Stabmittelfaser durch die 
jeweiligen Verschiebungen « in tangentialer Richtung 
und in Normalenrichtung w ausdrücken: 


1 IN 8 
at R, ) N as (3). 


Nimmt man wie üblich dehnungslose Verformung der 
Mittelfaser an, so ist weiter u"—=w und aus (3) folgt 
dann mit der Grundgleichung der Biegelehre 


8 
mRo 
EJ 


u"’+u=R, dd’ +cRo=f(P). . (4), 
0 


worin c die Integrationskonstante bezeichnet. Die 
vollständige Lösung von (4) ergibt die Tangential- 


. verschiebung (9): 


. 2 vo 
u(9) —sin 8 [ cos 9 f(d") a9 
Dr 
» . (5), 
— co8 »[ sin 9 f(9') 49’ 
Ds 


wenn d®, und 9, willkürliche Werte zur Festlegung von 
Integrationskonstanten bedeuten, 

Da u ebenfalls die Periode 2x besitzen muß, hat 
man die von m(d) abhängige Lösung (5) durch ge- 
eignete Wahl von m entsprechend zu bestimmen. Es 
soll also u(d) =u(9® +2) für jedes ® gelten und 
diese Periodizitätsforderung verlangt nach (4) und (5): 


9+2n Ey en 
3 [7 m r Y 
sind f co8 8 (* a + ne)an 
Ö ö 


+2n on R 
— cos # [ sin [m] a + m)anro 
0 


dv 


Die Koeffizienten von sin ® und cos 9 müssen einzeln 
verschwinden. 


Weil ferner der Krümmungsradius R, (#) wegen der 
Symmetrie der Ellipse die Periode x hat, verschwinden 
die Integrale 

+27 +27 
[ Recos®a®— | Rsind'dd =0 
» 02 


und als zweite und dritte „Randbedingung‘‘ des Eigen- 
wertproblemes (2) sind die Integralbedingungen 


27 d 
: [Rsins(/mz, a \a0=0 RT 
0 0 


en 02 
[Roco o(|ma,a0‘) ao NR ea 
0 0 R 


für m (®) gefunden. In (2)’’ und (2)’’” ist die Inte- 
gration über die Periode ohne Einschränkung der 
Allgemeinheit zwischen den Grenzen 0 und 2 aus- 
geführt und die Biegesteifigkeit #.J als konstant vor- 
ausgesetzt. 


4.Koeffizienten der Differentialgleichung. 
Die jetzt zu berechnenden Koeffizienten der Diffe- 
rentialgleichung (2) werden von der periodischen 
Funktion R,(®) und ihrem Kehrwert gebildet. Die 
Ellipsenkrümmung wird üblicherweise in Abhängig- 
keit von der Anomalie’i angegeben [3]: 
1 b — 3/2 
er RL) 
R, „(ilre cos? t) N 
worin & die durch b? — a = a?e? definierte zweite 
numerische Exzentrizität bedeutet. Mit atgt = btg ® 
läßt sich der Normalenwinkel gegen die + x-Achse 
einführen, und man erhält nach kurzer Rechnung 


&? cos? X 


Somit kann man in den Fourierreihen 


mn D} gm ©0852 md, ge D) Cam c08 2 md 
Ro m=0 7 m=0 


die Entwicklungskoeffizienten durch Potenzreihen 


von e darstellen, und man findet, daß cm und Cm 


jeweils proportional der Potenz e?m beginnen. Ins- 
besondere ist für den Grenzfall des Kreises (e=0) nur 


«=(,=1 von Null verschieden und für fast kreis- 


förmige Ellipsen ist e sehr klein, so daß es nahe liegt, 
nach Art der Störungsrechnung bei der Lösung des 
Eigenwertproblemes vorzugehen. Jedoch ist für 
praktisch in Frage kommende Achsenverhältnisse b/a 
etwa 1,5 e>1 und nicht hinreichend klein für eine 
ausschließliche Behandlung mit Hilfe der Störungs- 
rechnung. In der numerischen Rechnung sind daher 
die Fourierkoeffizienten 03m und 03m durch harmo- 
nische Analyse bequem zahlenmäßig bestimmt 
worden. 


5. Lösung des Eigenwertproblemes durch 
speziellen Reihenansatz 

Da die Differentialgleichung (2) von ungerader Ord- 
nung ist und überdies die Randbedingungen (2)” 
und (2)’’’ Integrale enthalten, sind die weiterreichen- 
den Methoden der Eigenwertberechnung bei dem vor- 
liegenden Problem nicht anwendbar. Daher werden 
die Eigenfunktionen in Fourierreihen entwickelt und 
der Eigenwert aus einer Abschnittsdeterminante des 
entstehenden unendlichen Gleichungssystems für die 
Fourierkoeffizienten genähert berechnet. — Die sämt- 
lichen Randbedingungen erfüllt der Ansatz 


m (d) = 3 (%n c0s2nd + Pansin2nd) 0), 


n=0 


denn die Periodizitätsbedingung (2)’ ist für jeden 


Fourieransatz erfüllt. Auch (2)” und (2) ist be- 
friedigt, da der Ansatz (6) verhindert, daß Quadrate 
über trigonometrische Funktionen in den Integranden 
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von (2)’’ und (2)’’’ vorkommen, und daher verschwin- 
den alle Integrale über die Periode 2 x einzeln. 

Von der Näherungsrechnung wird noch zusätzlich 
gefordert, daß sie stetig mit e den bekannten kleinsten 
Eigenwert für die Ausartung der Ellipse zum Kreis 
(e=0) richtig liefert. Für e=0 soll also aus den 
unendlichen Determinanten für die Fourierkoeffi- 
zienten A= — 3 folgen. 

Da die Differentialgleichung (2) linear ist, kann man 
die sich ergebenden Gleichungssysteme für die Fourier- 
koeffizienten &„ und fan getrennt betrachten. Durch 
Koeffizientenvergleich findet man mit dem Ansatz (6) 
zunächst für die ß,„ nach Einsetzen in (2): 


cos0 df,ß, +2c,ßı + 30, Ps -- -} 
+ 00829 (6, +34 +4(20, — Cy)]Pz 
+[66,+60+4(0,— C,)] By 
+94 +9Io +0, —Co)1Pe''} 
+ 084% {[l50,+ 150, +4(20,— 2C,)] Pa 
+[60 0, + 3005, +A[4C, — Ca)] Pa 
+[4505+ 45004420, —2Cy)]Ps' 30. 


Setzt man in der zugehörigen Koeffizientendeter- 
minante e=(, so verschwinden alle Elemente der 
ersten Zeile und alle Nichthauptdiagonalelemente. Es 
ergibt sich somit kein eindeutiger Eigenwert für den 
Grenzfall des Kreises, und diese Näherung wird daher 
als unbrauchbar ausgeschieden. 


Für die in bezug auf die Halbachsen’ der Ellipse 
symmetrischen Beulformen gewinnt man nach einem 
Koeffizientenvergleich das homogene Gleichungs- 
system, 


oo 
Don {n(4k®—1)[ogy on + gn2E Cgn+gzl 
n2=1 


+Ak [gg ont Conax t Contaz]}sin219 =0 
Br 7) 


für die Fourierköeffizienten An. In (7) sind &gm 
und C,m gleich Null zu setzen, wenn negative Indizes 
auftreten. Die unendliche Determinante der Koeffi- 
zienten &%n hat für &=0 nur Hauptdiagonalelemente 
und liefert, wie verlangt, dann auch den Eigenwert 
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)= — 3. Die Koeffizientendeterminante von (7) wird 
daher als brauchbare Näherung zur Berechnung des 
kleinsten Eigenwertes angesehen. Für verschiedene 
Werte des Achsenverhältnisses b/a wurden die Fourier- 
koeffizienten (3m und C,m zahlenmäßig bestimmt und 
dann ein Näherungswert A durch die Nullstellen der 
Abschnittsdeterminante von (7) berechnet. Bei dieser 
Determinantenauswertung konnte man sich für die 
vorliegende technische Fragestellung auf zwei- oder 
dreireihige Abschnittsdeterminanten mit hinreichend 
genauen Ergebnissen beschränken. 


6. Ergebnisse 
Die numerische Rechnung lieferte die Zahlenwerte: 


.=fire 1.001.110 4°1,20° 1.38.) 1.00.0708 


Sa? 
Mi =3.00 2.70 22,40 —2.00.-1,0 
ET 3,00 —2,70 —2,40 —2,00 —1,70 —1,45 
N er 5 
Damit ist die kritische Belastung Z; = — Fr: 


EJ 1 

: ® BR, 
geben. Über die Güte und Brauchbarkeit des Reihen- 
ansatzes für m(d) lassen sich keine allgemeinen 
Aussagen machen. So kann nur die Zahlenrechnung 
wenigstens für kleine e geprobt werden. Mit Hilfe der 


abhängig vom Achsenverhältnis a ange- 
a 


Störungsrechnung findet man nämlich = — 3 -+ 5 & 


..., was für sehr kleine e Übereinstimmung mit den 
auf anderem Wege berechneten Zahlen ergibt. 
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BUCHBESPRECHUNGEN 


Prof. Dr. L. Bieberbach, Einführung in die 
Analytische Geometrie. Vierte, neubearbeitete 
Auflage. 168 S. mit 43 Abb. Bielefeld 1950. Verlag 
für Wissenschaft und Fachbuch GmbH. Preis: 
kart. 8,90 DM. 

Das Buch ist eine Neubearbeitung der bewährten, 
früher im Teubnerverlag erschienenen ‚‚Analytischen 
Geometrie“ des gleichen Verfassers; es soll deshalb 
hauptsächlich über die vorgenommenen Verände- 
rungen berichtet werden. Das Charakteristische des 
Buches ist durchaus gewahrt geblieben: Vf. gibt auf 
knappem Raum eine Einführung in die Analytische 
Geometrie und stellt gleichzeitig die methodischen 
Grundlagen, Vektorrechnung und Matrizenkalkül, 
dar. Eine scharfe Trennung zwischen beiden Gebieten 
. ist jedoch nicht gemacht, vielmehr werden jeweils 
nach Bereitstellung der nötigen Begriffe und Formeln 
(für den n-dimensionalen Raum) wieder Beispiele und 
Anwendungen geometrischer Art gebracht und so die 
analytische Geometrie gewissermaßen am Rande auf- 
gebaut. Gegenüber früheren Auflagen fällt vor allem 
auf, daß von der dort angedeuteten Möglichkeit, auf 
den Begriff des Ortsvektors zu verzichten, konsequent 
Gebrauch gemacht wird. Punkte und Vektoren im R, 
werden gemeinsam durch eine Folge von n+1 Zahlen 
dargestellt, deren erste n mit den üblichen Koor- 
dinaten übereinstimmen, während die (n-+1)-te für 


Punkte gleich 1, für Vektoren gleich O ist. Sämtliche 
Rechnungen werden mit derartigen Zahlenfolgen aus- 
geführt, in denen Punkte (an Stelle der Ortsvektoren) 
und Vektoren gleichberechtigt auftreten. Diese me- 
thodische Einfachheit wird allerdings auf Kosten eines 
zunächst ungewohnten Schriftbildes erzielt. Neben 
verschiedenen Änderungen und Umstellungen, z. B. 
einer ausführlicheren Darstellung der Determinanten 
ist die Neuauflage um zwei Abschnitte bereichert wor- 
den, die dem Begriff des Flächeninhalts und der 
Kreisgeometrie gewidmet sind. Dank seiner besonderen 
Anlage gibt das Buch dem aufmerksamen Leser nicht 
nur einen ausgezeichneten ersten Einblick in die 
analytische Geometrie und die lineare Algebra, son- 
dern darüber hinaus ein schönes Beispiel mathema- 
tischer Denkweise. 


Dresden. G. Opitz. 


W. Lietzmann, (Prof. a. d. Univ. Göttingen), Me - 
thodik des mathematischen Unter- 
riehts (der Lehrstoff), 2088. Heidelberg 1951. 
Quelle u. Meyer. Preis geb. 13,80 DM. 

Es handelt sich um eine den Zeitverhältnissen an- 
gepaßte gekürzte Neubearbeitung der vom Vf, früher 
veröffentlichten und weit verbreiteten Methodik. Der 
Wert des Buches liegt wesentlich inseiner vorzüglichen 
Eignung als Nachschlagewerk. Es dürfte außer 
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Lietzmann kaum jemand eine so umfassende 
Kenntnis der einschlägigen Literatur und aller Ent- 
wicklungsphasen des mathematischen Unterrichts be- 
sitzen. Dadurch kommt allerdings der Standpunkt 
des V£. nicht hinreichend klar zum Ausdruck. Das ist 
vor allem deshalb zu beklagen, weilso die grundsätz- 
liche Stellungnahme zu denjenigen Ansprüchen an 
die Schule unbestimmt bleibt, diesich ausder modernen 
Entwicklung der Wissenschaft und des Lebens unab- 
weisbar ergeben. In dieser Hinsicht bedarf die junge 
Generation, der das Buch gewidmetist, deserfahrenen 
methodischen Rates wohl doch notwendig. 

Zur Erläuterung seien hier die Ausführungen zur 
Wahrscheinlichkeitsrechnung wiedergegeben. Zu- 
nächst wird die Einführung in die „traditionelle 
Wahrscheinlichkeitslehre‘‘ empfohlen, und später 
heißt es: „Es ist zu wünschen, daß man mit der 
neueren Grundlegung Unterrichtsversuche macht. 
Aber auch wer an der alten Einführung festhält, sollte 
der neueren Auffassung Rechnung tragen‘. Über das 
Wie der praktischen Verwirklichung dieser Ratschläge 
wird nichts gesagt, und das ist ohne Erfahrung nicht 
zuleisten. Die zitierten Empfehlungen führen, wenn 
man sie einmal konkret durchdenkt, entweder zur Be- 
griffsverwirrung oder in eine Diskussion der Grund- 
lagenkrise hinein, die wissenschaftlich exakt zu führen 
— und jedes oberflächliche Hinweggleiten verbietet 
dem Mathematiker der Geist seiner Wissenschaft —, 
über den Rahmen der Schule hinausgeht. Im übrigen 
kann auch die Schule — und diese Forderung wäre 
angesichts der heutigen Bedeutung der Wahrschein- 
lichkeitsrechnung an sie zu stellen — auf die klassische 
Lehre verzichten und dem Unterricht z. B. die Auf- 
fassung von v. Miseszugrundelegen. f 

Dresden. Draeger. 

William T. Sawyer, Experimental Investi- 
gation of a Stationary Cascade of Aero- 
dynamic Profiles. (Mitt. Institut f. Aerodynamik 
an der E. T.H. Zürich, herausgegeben von Prof. 
Dr. J. Ackeret, Nr.17) 788. Zürich 1949. Verlag 
Leemann. Preis 12,— Fr. 

In dieser beachtlichen Forschungsarbeit wird die 
Ackeretsche Theorie für den Entwurf dichtstehender 
Schaufelgitter experimentell geprüft. Um die Strö- 
mung durch ein solches Sohrutelgitten mit vor- 
geschriebenen Wirkungen zu berechnen, nimmt 
Ackeret periodische kontinuierliche Verteilungen 
von Wirbeln und Quellen an, dieam Ort der Schaufeln 
stark konzentriert gedacht werden. Das Schaufel- 
profil, das in dem jetzt untersuchten einzigen Fall 
herauskommt, weist eine gewünschte Ähnlichkeit mit 
dem N.A.C.A.-Profil 65—010 auf. Die Luftkräfte bei 
Geschwindigkeiten, die weit unter der Schallgeschwin- 
digkeit liegen, werden mit einer Spezialwaage ge- 
messen, die eigens für Flügelgitter entwickeltist. Um 
den Einfluß der Reynoldsschen Zahl und den der 
Turbulenz des Luftstromes zu übersehen, wird, ein 
breit angelegtes Versuchsprogramm durchgeführt, 
wobei sich zeigt, daß die Intensität der Turbulenz an 
Einfluß weit den der Skala übertrifft. Diese beiden 
Eigenschaften werden im Anschluß an die von 
G. I. Taylor inaugurierten Untersuchungen zur 
Turbulenz von Windströmen bestimmt. Wenn die 


a En; N 


Z.angew. Math. Mech. 
Bd. 31 Nr.10 Okt.1951 


Reynoldssche Zahl nicht zu klein ist, so ergibt 
sich eine vorzügliche Übereinstimmung der Ackeret- 
schen Theorie mit diesen Experimenten. 

Zweifellos wird diese sorgfältige und besonnene 
Arbeit sowohl dem Strömungsforscher als auch dem 
Konstrukteur von Strömungsmaschinen wertvolle 
Einsichten vermitteln. 


Göttingen. W. Tollmien. 


Dr. Horst Teichmann (Universität Würzburg), Ein- 
führung in die Quantenphysik. (Mathematisch- 
Physikalische Bibliothek, Reihe II, Bd. 13) 2. Aufl., 
104 S. mit 44 Abb. Leipzig 1950. B. G. Teubner 
Verlagsgesellschaft. Preis: kart. 4,20 DM. 

Das Büchlein stellt die experimentellen Grundlagen 
der Quantenhypothese dar und bringt die wichtigsten 
Anwendungen dieser Theorie. Sowohl die Physik der 
Atomhülle als auch die des Atomkernes wird be- 
handelt. Das Buch ist für eine breite Schicht natur- 
wissenschaftlich Interessierter gedacht. Dement- 
sprechend ist die Darstellung leicht verständlich, 
trotzdem aber frei von dem Bemühen um unwissen- 
schaftliche Popularisierung des behandelten Stoffes. 


A. Recknagel. 


P. Crantz. Sphärische Trigonometrie. Bearb. 
von Dr. M. Hauptmann. (Mathematisch-physika- 
lische Bibliothek, Reihe II, 7.) 4. Aufl., 1128. m. 
76 Abb. und 1 Nomogramm. Leipzig 1950. B. G. 
Teubners Verlagsgesellschaft. Preis kart. 4,20 DM. 

In ausführlicher und anschaulicher Weise, die für 
den Selbstunterricht geeignet ist, wird, von gering- 
fügigen Erweiterungen abgesehen, der'in der Schule 
zu behandelnde mathematische Stoff nebst seinen 
Anwendungen auf Erd- und Himmelskunde darge- 
stellt. Auch der Zusammenhang zwischen sphärischer 
und ebener Trigonometrie und die regulären Polyeder 
werden abgeleitet. Eine große Zahl von Übungs- 
aufgaben ermöglicht dem Leser ein gründlicheres 
Vertrautwerden mit dem Stoff. Nur ist leider wegen 
ungenauer Wiedergabe das Nomogramm (Fig. 46, 47) 
nicht benutzbar. 

Dresden. 


Dr.-Ing. Hans Mahl und Erich Gölz. Elektronen- 
mikroskopie. (Die kleinen Studienbücher, Band58) 
183 S. 88 Abb. und 39 Aufnahmen. Leipzig 1951. 
VEB Bibliographisches Institut. Preis kart. 5,80 DM. 

Inhalt: Grundlagen der Elektronenmikroskopie. 
Aufbau und Eigenschaften der verschiedenen Elek- 
tronenmikroskope. Konstruktionsteile und Zusatz- 
einrichtungen. „Arbeitsmethoden und praktischer 
Einsatz des Übermikroskops. 

In dem Rahmen, der durch diese Stichworte ge- 
kennzeichnet ist, bringen die Verfasser einen voll- 
ständigen Abriß der Elektronenmikroskopie. Dabei 
wenden sie sich vor allen Dingen an Nichtphysiker, 
die sich in die Grundlagen der Elektronenmikroskopie 
einarbeiten müssen, ohne daß sie die umfangreichen, 
für Fachleute gedachten Darstellungen studieren 
könnten. Das Buch löst seine Aufgabe mit Geschick 


Draeger. 


und wird sich viele Freunde erwerben. Der niedrige 


Preis ist besonders hervorzuheben. 


Dresden. Recknagel. 


NACHRICHTEN 


Halle (Saale): Am 7. Sept. 1951 starb der o. Prof. 
der angewandten Mathematik an der Universität Halle 
Dr. Harry Schmidt. 


Mainz:Dr.HeinrichHenckyverunglückte 
am 6. Juli 1951 tödlich bei Ausübung des Bergsportes. 
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